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^ Abstract 

^ I present a pedagogical survey of a variety of quantum phases of the Hubbard model. The 

C^l honeycomb lattice model has a conformal field theory connecting the semi-metal to the insulator 

with Neel order. States with fractionalized excitations are linked to the deconfined phases of 
gauge theories. I also consider the confining phases of such gauge theories, and show how Berry 
^ phases of monopoles induce valence bond solid order. The triangular lattice model can display 

a metal-insulator transition from a Fermi liquid to a deconfined spin liquid, and I describe the 

^ theory of this transition. The bilayer triangular lattice is used to illustrate another compressible 

metallic phase, the 'fractionalized Fermi liquid'. I make numerous connections of these phases and 
critical points to the AdS/CFT correspondence. In particular, I argue that two recent holographic 

^ constructions connect respectively to the Fermi liquid and fractionalized Fermi liquid phases. 
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I. INTRODUCTION 



The Hubbard model is the simplest of a class of models describing electrons moving on 
a lattice with repulsive electron-electron interactions. Despite its apparent simplicity, it has 
become clear in the past two decades that it can display a very rich phase diagram, with a 
plethora of interesting phases. The most common phase is, of course, the Fermi liquid (FL), 
which is adiabatically connected to the metallic phase of non-interacting electrons. However, 
electron-electron interactions can break one or more symmetries of the Hamiltonian, and this 
leads to phases such as antiferromagnets, charge or spin density waves, or superconductors. 
Also of great interest are quantum phases which do not break any symmetries, but are 
nevertheless qualitatively distinct from the non-interacting electron states: such states are 
characterized by emergent gauge excitations, fractionalization of quasiparticle excitations, 
and non-trivial ground state degeneracies which depend upon the global topology of the 
lattice — it is often stated that such states have 'topological' order. Finally, there are inter- 
esting quantum phase transitions between such phases, and such quantum critical points 
are often described by strongly-coupled quantum field theories. In some cases, the quantum 
critical points can broaden into gapless quantum critical phases. 

This article will present a pedagogical review of a small sample of this landscape of phases 
and critical points. My aim is to describe the appearance of a variety of non-trivial phases in 
the simplest possible context. For the honeycomb lattice with a density of one electron per 
site, such phases naturally have low energy excitations which have a relativistic form at low 
energies. Consequently, in the vicinity of quantum phase transitions, such phases and their 
critical points are amenable to a description by relativistic quantum field theories. In some 
cases, the critical points are also conformally invariant, and so are described by conformal 
field theories (CFTs). In these cases, the AdS/CFT correspondence can be directly applied, 
and I will describe the insights that have been gained from such an approach. 

However, once we move away from commensurate electron densities, the quantum phases 
and critical points of electron lattice models rarely have any relativistic invariance in their 
low energy theory. I will describe here the simplest examples of 'topologically ordered' 
phases at generic electron densities. It is important to note that the lack of relativistic 
invariance does not rule out application of the AdS/CFT correspondence. We can begin 
from a relativistically invariant gravity dual theory and dope it with charge carriers by 
turning on a chemical potential: then even the gravity theory is not relativistically invariant 
at low energies, and we can hope to match its low energy physics to a condensed matter 
system. There has been a large effort to apply the AdS/CFT correspondence along this 
direction in the past few years. I have discussed some of this work in another recent review 
articl^, which should be viewed as a companion to the present article. I will also discuss 
the generic density phases here, with an emphasis on the general low energy structure of the 
'fractionalized Fermi liquid' (FL*) phas^^, which I believe is closely related to the generic 
density phases that have appeared using the AdS/CFT correspondence. Further details on 
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the connection of these phases of Hubbard-hke models and the AdS/CFT correspondence 
appear in a recent papei'^. 

I begin by introducing the Hubbard model. It is defined by the Hamiltonian 



^ = -Yl + 5Z ~^ ^^^^ + ^^^^ + (^*t - \^ (^4 - \^ 



Here Qq,, a =t, i are annihilation operators on the site i of a regular lattice, and tij is 
a Hermitian, short-range matrix containing the 'hopping matrix elements' which move the 
electrons between different lattice sites. The density of electrons is controlled by the chemical 
potential ^ which couples to the total electron density, with 

n,i^ = cJ^Qt , = cl^di. (1.2) 

The electrons repel each other with an on-site interaction Uf, in most cases we will take 
Ui = U site-independent, but it will also be useful later to allow for a site-dependent Ui. For 
completeness, we also note the algebra of the fermion operators: 

CiaCjp + CjpCia =0. (1-3) 



The equations (1.1), (1.2), and (1.3) constitute a self-contained and complete mathematical 
statement of the problem of the landscape of the Hubbard model. It is remarkable that a 
problem that is so simple to state has such a rich phase structure as a function of the lattice 
choice, the fermion density, and the spatial forms of tij and f/j. 

Sections [IT1,|IV[ and |V] will deal exclusively with the honeycomb lattice at a density of one 
electron per site ("half-filling"), so that (?T-jt) = (^4) = 1/2. The emphasis on the honeycomb 
lattice is not motivated by its particular physical importance (although, it is the lattice of 
graphene), but by its simplicity as a context for introducing various technical methods, 
quantum phases and critical points. In Section [TT| we will consider the semi-metal and the 
insulating antiferromagnet, and show that a phase transition between them is described by a 



relativistic field theory, which is a version of the Gross- Neveu- Yukawa model; Section HI will 
use this field theory to present a general discussion of the physics at non-zero temperatures 
in the vicinity of a quantum critical point. We will focus on the transport of conserved 



charges, and describe insights gained from the AdS/CFT correspondence. Section IV will 
consider the problem of restoring the spin rotation symmetry from the antiferromagnet, 
while remaining in an insulating phase: this will lead to a description in terms of a U(l) 
gauge theory, and the appearance of an insulating phase with valence bond solid (VBS) 
order. Finally, Section |V] will combine all the phases of the half-filled honeycomb lattice 
discussed so far in a single phase diagram: this will require introduction of a SU(2) gauge 
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FIG. 1: The honeycomb lattice with its A (green) and B (red) sublattices 



theory. We will find an interesting multi-critical point in Section |V| which has many features 
in common with the supersymmetric CFTs studied using the AdS/CFT correspondence. 



Sections IV and Mean be skipped in a first reading, without significant loss of continuity. 



We will move away from half-filling in Sections VI and VII, where we will turn our 



attention to metallic phases with Fermi surfaces. Section |VI] considers the Hubbard model 
on the triangular lattice, and describes a phase diagram which includes Fermi liquid (FL) 



and spin liquid phases. Section |VII| extends our discussion to the the Hubbard model on 
a bilayer triangular lattice, which has been realized in experiment^ on ^He. Here we will 
present a gauge theory of another metallic phase, the fractionalized Fermi liquid (FL*). 
We will also discuss the connections to compressible metallic phases obtained using the 
AdS/CFT correspondence, complementing the recent discussion in Ref. O 



II. SEMI-METAL AND ANTIFERROMAGNETISM ON THE HONEYCOMB 
LATTICE 

A. Preliminaries 



We will consider the Hubbard model (1.1) with the sites i on locations on the hon- 



eycomb lattice shown in Fig. [T} Here, we set up some notation allowing us to analyze the 
geometry of this lattice. 

We work with a lattice with unit nearest neighbor spacing. We define unit length vectors 
which connect nearest-neighbor sites 



ei = (l,0) , e2 = (-l/2,v^/2) , eg = (-1/2, -^3/2). 



(2.1) 



Note that ■ e,- = —1/2 for i ^ j, and 61 + 62 + 63 = 0. The lattice can be divided into the 
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FIG. 2: The first Brillouin zone of the honeycomb lattice. 



A and B sublattices, as shown in Fig. [T] We take the origin of co-ordinates of the lattice at 
the center of an empty hexagon. The A sublattice sites closest to the origin are at ei, 62, 
and 63, while the B sublattice sites closest to the origin are at — Ci, —62, and —63. 

The unit cell of the hexagonal lattice contains 2 sites, one each from the A and B sub- 
lattices. These unit cells form a triangular Bravais lattice consisting of the centers of the 
hexagons. The triangular lattice points closest to the origin are ±(ei — 62), ±(e2 — 63), and 
±(63 — Ci). The reciprocal lattice is a set of wavevectors G such that G ■ r = 27tx integer, 
where r is the center of any hexagon of the honeycomb lattice. The reciprocal lattice is also 
a triangular lattice, and it consists of the points '^^UiGi, where are integers and 



An 
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G, 



An 



-63- 



(2.2) 



The unit cell of the reciprocal lattice is called the first Brillouin zone. This is a hexagon 
whose vertices are given by 



Qi = liG2-Gs) , Q2 = liGs~G,) 



Q, = ^(G,-G2) 



(2.3) 



and —Qi, —Q2, and —Q3] see Fig. |2j Integrals and sums over momentum space will 
implicitly extend only over the first Brillouin zone. This is the 'ultraviolet cutoff' imposed 
by the underlying lattice structure. 

We define the Fourier transform of the electrons on the A sublattice by 



(2.4) 



where M is the number of sites on one sublattice; similarly for csa- Note that CAa{k + G) = 
CAa{k): consequently, sums over momentum have to be restricted to the first Brillouin zone 



to avoid double counting. Thus the inverse of Eq. (2.4) sums over k in the first Brillouin 
zone. 
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FIG. 3: The lower band of the dispersion in Eq. (2.7) for n = 



B. Semi-metal 



We begin with free electrons in the honeycomb lattice, U = 0, with only nearest-neighbor 
electron hopping tij = t. Using Eq. (2.4), we can write the hopping Hamiltonian as 



Ho = -tJ2{e''''"' + e*''"' + e*'') claik)cBa{k) + H.c. 

k 



(2.5) 

We introduce Pauli matrices r" (a = x, y, z) which act on the A, B sublattice space; then 
this Hamiltonian can be written as 

Hq = ''^c^{k) — /X — t^cos(A; ■ ei) + cos(fc ■ 62) + cos(fc ■ 63) jr^ 

k 

+ t(^sm{k ■ ei) + sin(fc • 62) + sin(fe ■ 63) jr^ c(fc), (2.6) 

where the sublattice and spin indices on the electrons are now implicit: the c{k) are 4- 
component fermion operators. 

The energy eigenvalues are easily determined to be 



(2.7) 



and these are plotted in Fig. [3j At half-filling, exactly half the states should be occupied in 



the ground state, and for the spectrum in Eq. (2.7) this is achieved at /i = 0. 



A crucial feature of any metallic state is the Fermi surface: this is boundary between the 
occupied and empty states in momentum space. In two spatial dimensions, this boundary 



is generically a line in momentum space, and this is the case for the dispersion in Eq. (2.7) 
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for n ^ 0. However, for the /i = 0, the honeycomb lattice has the special property that the 
occupied and empty states meet only at a discrete set of single points in momentum space: 
this should be clear from the dispersion plotted in Fig. |3} Only 2 of these points are distinct, 
in that they are not separated by a reciprocal lattice vector G. So the half-filled honeycomb 
lattice has 2 'Fermi points', and realizes a 'semi-metal' phase. The low energy excitations 
of the semi- metal consist of particles and holes across the Fermi point, and these have a 
lower density of states than in a metallic phase with a Fermi line. We also note that the 
Fermi-point structure is protected by a sublattice exchange symmetry: it is not special to 
the nearest-neighbor hopping model, and it also survives the inclusion of electron-electron 
interactions. 

We obtain a very useful, and universal, theory for the low energy excitations of the semi- 



metal by expanding (2.6) in the vicinity of the Fermi points. The distinct Fermi points are 
present at Qi and —Qi] all other Fermi points are separated from these two points by a 
reciprocal lattice vector G. So we define continuum Fermi field which reside in 'valleys' in 
the vicinity of these points by 



CAla{k) = VAcAa{Ql + k) 

CA2a{k) = VIcAai-Ql + k) 

Ceiaik) = VAcBaiQi + k) 

CB2aik) = ^CBci-Ql + k), (2.8) 

where A is the total area of the honeycomb lattice, and the momentum k is small. The field 
C is a 8-component continuum canonical Fermi field: the components correspond to spin 
(t, i), sublattice {A, B), and valley (1, 2) indices. We will also use Pauli matrices which act 
on the spin (a"), sublattice (r"), and valley {p") space. 



Inserting Eq. (2.8) into Eq. (2.6), we obtain the continuum Hamiltonian 



Ho = J ^C\k)(vTyK + vT^p'ky)Cik), (2.9) 



where v = 3t/2. From now on we rescale time to set v = 1. Diagonalizing Eq. (2.9), we 
obtain the relativistic spectrum 

± Jkl + kl (2.10) 



which corresponds to the values of Eq. (2.7) near the Fermi points. 

The relativistic structure of Hq can be made explicit by rewriting it as the Lagrangian 
of massless Dirac fermions. Define C = C'^ p^r^ . Then we can write the Euclidean time (r) 
Lagrangian density of the semi-metal phase as 

£o = C{drio + 5x71 + -9,72) C (2.11) 
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where oj is the frequency associated with imaginary time, and the Dirac 7 matrices are 

7o = -pV^ 1i = P't^ 72 = -r^. (2.12) 

In addition to relativistic invariance, this form makes it clear the free-fermion Lagrangian 
has a large group of 'flavor' symmetries that acts on the 8x8 fermion space and commute 
with the 7 matrices. Most of these symmetries are not obeyed by higher-order gradients, or 
by fermion interaction terms which descend from the Hubbard model. 

Let us now turn on a small repulsion, U , between the fermions in the semi-metal. Because 
of the point-like nature of the Fermi surface, it is easier to determine the consequences 
of this interaction here than in a metallic phase with a Fermi line of gapless excitations. 
We can use traditional renormalization group (RG) methods to conclude that a weak U is 
irrelevant in the infrared: the computation is left as an exercise below. Consequently, the 
semi-metal state is a stable phase which is present over a finite range of parameters. 

Exercise: Observe that £0 is invariant under the scaling transformation x' = xe~^ 
and r' = re~^. Write the Hubbard interaction U in terms of the Dirac fermions, and show 
that it has the tree-level scaling transformation U' = Ue~^. So argue that all short-range 
interactions are irrelevant in the Dirac semi-metal phase. 



C. Antiferromagnet 

Although a small U is irrelevant, new phases can and do appear at large U. To see 



this, let us return to the lattice Hubbard model in Eq. (1.1), and consider the limit of large 
Ui = U. We will assume /i = and half-filling in the remainder of this section. 

At t/ = 00, the eigenstates are simple products over the states on each site. Each site 
has 4 states: 

|0) , 4^10) , 4J0) , 4^4J0), (2.13) 

where |0) is the empty state. The energies of these states are U/A, — f//4, — f//4, and 11/4 
respectively. Thus the ground state on each site is doubly-degenerate, corresponding to the 
spin-up and spin-down states of a single electron. The lattice model has a degeneracy of 2^-^, 
and so a non-zero entropy density (recall that JV is the number of sites on one sublattice). 

Any small perturbation away from the U = 00 limit is likely to lift this exponential 
large degeneracy. So we need to account for the electron hopping t. At first order, electron 
hopping moves an electron from one singly-occupied site to another, yielding a final state 
with one empty and one doubly occupied site. This final state has an energy U higher than 
the initial state, and so is not part of the low energy manifold. So by the rules of degenerate 
perturbation theory, there is no correction to the energy of all the 2^-^ ground states at first 
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order in t. 

At second order in t, we have to use the effective Hamiltonian method. This performs 
a canonical transformation to ehminate the couphngs from the ground states to all the 
states excited by energy U, while obtaining a modified Hamiltonian which acts on the 2^-^ 
ground states. This method is described in text books on quantum mechanics, and we leave 
its application here as an exercise. The resulting effective Hamiltonian is the Heisenberg 
ant if erromagnet : 

HJ = J2■J^^StSJ , ■J^, = -^, (2.14) 

i<j 

where Jij is the exchange interaction and S!^ are the spin operators on site i 

St = \cl<pc.p. (2.15) 

Note that these spin operators preserve the electron occupation number on every site, and 
so act within the subspace of the 2"^ low energy states. The Hamiltonian Hj lifts the 
macroscopic degeneracy, and the entropy density of the new ground state will be zero. 



Exercise: Use the effective Hamiltonian method described in Ref. [7] to obtain Eq. (2.14). 
At second order in f/, it is sufficient to consider the 2-site Hubbard model. This has a total 
of 16 states, and 4 ground states ai U = oo. Derive the effective Hamiltonian which acts 
on these 4 states at order t^. 



Although we cannot compute the exact ground state of Hj on the honeycomb lattice 
with nearest-neighbor exchange, numerical studied leave little doubt to its basic structure. 
The ground state is adiabatically connected to that obtained by treating the S'f as classical 
vectors in spin space: it has antiferromagnetic (or Neel) order which breaks the global SU(2) 
spin rotation symmetry, by a spontaneous polarization of the spins on opposite orientations 
on the two sublattices 

T^,{St) = N\ (2.16) 

where r]i = 1 {rji = — 1) on sublattice A (B), and A^" is the vector Neel order parameter; see 
Fig. [4j Classically this state minimizes the exchange coupling in Eq. (2.14) because Jij > 0. 



Quantum fluctuations for spin 5 = 1/2 reduce the spontaneous moment from its classical 
value, but a non-zero moment remains on the honeycomb lattice. 

What is the electronic excitation spectrum in the ant if erromagnet? To determine this, it 
is useful to write the Neel order parameter in terms of the continuum Dirac flelds introduced 
in Section Til B[ We observe 

E "1^^^ = E {^Aa<pCA, - C^Ba<pCB,) = [ '^.C^^o'^C (2.17) 
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FIG. 4: The large U state with antiferromagnetic (Neel) order. 



Thus the Neel order parameter N°- is given by the fermion bihnear 

= = (CpV"C) , (2.18) 

and the vacuum expectation value (VEV) is non-zero in the antiferromagnet. We can expect 
that electron-electron interactions will induce a coupling between the fermion excitations and 
this VEV in the low energy Hamiltonian for the Neel phase. Choosing Neel ordering in the 
z direction 

iV^ = NoSa,, (2.19) 



we anticipate that Hq in Eq. (2.9) is modified in the Neel phase to 

Hn = J 0Ct(fc)(rJ^A;, + rVfc,-AiVor%^)c(fc), (2.20) 

where A is a coupling determined by the electron interactions, and we have assumed Neel 
order polarized in the z direction. This effective Hamiltonian will be explicitly derived in the 
next subsection. We can now easily diagonalize to deduce that the electronic excitations 
have energy 

± ^kl + kl + Xmi (2.21) 

This is the spectrum of massive Dirac fermions. So the Fermi point has disappeared, and an 
energy gap has opened in the fermion excitation spectrum. In condensed matter language, 
the phase with antiferromagnetic order is an insulator, and not a semi-metal: transmission 
of electronic charge will require creation of gapped particle and hole excitations. 
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D. Quantum phase transition 



We have now described a semi-metal phase for small f/, and an antiferromagnetic insulator 
for large U. Both are robust phases, whose existence has been reliably established. We now 
consider connecting these two phases at intermediate values of U. This is a complex subject, 
and careful numerical studies are only just emerging for the model with nearest-neighbor 
hopping. It is already clear, however, that by varying the form of the microscopic coupling 
we can obtain a rich variety of intermediate phase^^^. In the present subsection we consider 
the simplest possibility: there are no new intermediate phases, and only a direct quantum 
phase transition between the semi-metal and the antiferromagnetic insulatoi'^^'^. 

We can derive the field theory for this direct transition either by symmetry considerations, 
or by an explicit derivation from the Hubbard model. Let us initially follow the second route. 



We start with the Hubbard Hamiltonian in Eq. (1.1), use the operator identity (valid on 
each site i): 



U 



2U 



u 



(2.22) 



Then, in the fermion coherent state path integral for the Hubbard model, we apply a 
'Hubbard-Stratonovich' transformation to the interaction term; this amounts to using the 
identity 




drS: 



a2 



:X. 



a2 



ux^si 



(2.23) 



The fermion path integral is now a bilinear in the fermions, and we can, at least formally, 
integrate out the fermions in the form of a functional determinant. We imagine doing this 
temporarily, and then look for the saddle point of the resulting effective action for the X^. 
At the saddle-point we find that the lowest energy is achieved when the vector has opposite 
orientations on the A and B sublattices. Anticipating this, we look for a continuum limit in 
terms of a field y?" where 

Xt = r^,^'^ (2.24) 



Using Eq. (2.17), the continuum limit of the coupling between the field if"" and the fermions 



in Eq. (2.23) is given by 



Xi c]^(T^aC' 



I la a 



(2.25) 



From this it is clear that ip'^ is a dynamical quantum field which represents the fluctuations 
of the local Neel order, and 

oc N". (2.26) 
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Now we can take the continuum limit of all the terms in the coherent state path integral 
for the lattice Hubbard model and obtain the following continuum Lagrangian density 

£ = C^,d,C + \ [{d.v'^f + s^'^'] + ^ (^'^2) ' - X^'^Cp^a'^C (2.27) 

This is a relativistic quantum field theory for the 8-component fermion field C and the 
3-component real scalar related to the Gross- Neveu- Yukawa model. We have included 
gradient terms and quartic in the Lagrangian for <^°: these are not present in the derivation 
outlined above from the lattice Hubbard model, but are clearly induced by higher energy 
fermions are integrated out. The Lagrangian includes various phenomenological couplings 
constants (s, u, A); as these constants are varied, C can describe both the semi-metal and 
insulating antiferromagnet phases, and also the quantum critical point between them. 

Note that the matrix p^a"^ commutes with all the 7^; hence p^cr*^ is a matrix in "flavor" 
space. So if we consider C as 2-component Dirac fermions, then these Dirac fermions carry 
an additional 4-component flavor index. 

The semi-metal phase is the one where (p"" has vanishing VEV. In mean-field theory, this 
appears for s > 0. The if"" excitations are then massive, and these constitute a triplet of 
gapped 'spin-excitons' associated with fluctuations of the local antiferromagnetic order. The 
Dirac fermions are massless, and represent the Fermi point excitations of the semi-metal. 

The Neel phase has a non-zero VEV, 7^ 0, and appears in mean-field theory for 
s < 0. Here the Dirac fermions acquire a gap, indicating that the Fermi point has vanished, 
and we are now in an insulating phase. The fluctuations of (p are a doublet of Goldstone 
modes ('spin waves') and a longitudinal massive Higgs boson. 

Finally, we are ready to address the quantum critical point between these phases. In 
mean- field theory, this transition occurs at s = 0. As is customary in condensed matter 
physics, it is useful to carry out an RG analysis near this point. The tree- level analysis is 
carried out in the following exercise. 

Exercise: Perform a tree-level RG transformation on C The quadratic gradient 
terms are invariant under C = Ce^ and ip' = Lpe^l"^ . Show that this leads to s' = se^^. 
Thus s is a relevant perturbation which drives the system into either the semi-metal or 
antiferromagnetic insulator. The quantum critical point is reached by tuning s to its critical 
value (= at tree level). Show that the couplings u and A are both relevant perturbations 
at this critical point. Thus, while interactions are irrelevant in the Dirac semi- metal (and 
in the insulator), they are strongly relevant at the quantum-critical point. 

Further study of this quantum critical point requires a RG analysis which goes beyond 
tree- level. Such an analysis can be controlled in an expansion in (where N is the 
number of fermion fiavors) or (3 — d) (where d is the spatial dimensionality. For reviews 
see Ref. [18] or Chapter 17 of Ref. [191 The main conclusion of such analyses is that there is 
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FIG. 5: The honeycomb lattice with a vacancy. 



an RG fixed point at which the y?"^ is the only relevant perturbation. Non-linearities such 
as A and u all reach stable fixed point values of order unity. This non-trivial fixed point 
implies that the physics of the quantum critical point is highly non-trivial and strongly 
coupled. The RG fixed point is scale- and relativistic- invariant, and this implies that it is 
also conformally invariant. Thus the quantum critical point is described by a CFT in 2-1-1 
spacetime dimensions: a CFT3. 

We will not describe the critical theory in any detail here. However, we will note some 
important characteristics of correlation functions at the quantum critical point. The electron 
Green's function has the following structure 



where rj > is the anomalous dimension of the fermion. This leads to a fermion spectral 
density which has no quasiparticle pole: thus the quantum critical point has no well-defined 
quasiparticle excitations. This distinguishes it from both the semi-metal and insulating 
antiferromagnetic phases that fiank it on either side: both had excitations with infinitely- 
sharp quasiparticle peaks. Similar anomalous dimensions appear in the correlations of the 
bosonic order parameter ip"-. 

E. Quantum impurities 

We briefiy note the physics of quantum impurities in the honeycomb lattice, which were 
discussed more completely in the companion revieW^. The translational invariance of the 
honeycomb lattice will be broken only in this subsection. 

Imagine removing a single atom from the honeycomb lattice, as shown in Fig.|5]At U = 0, 
the electronic spectrum of such a lattice was described in Ref. |20l They showed that there 




iuj + kxT^ + kyT^'p^ 



(2.28) 
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was a quasi-localized state in the vicinity of the impurity exactly at zero energy, which 
decayed as 1/r at a distance r from the impurity. Upon including U, we expect this state 
to capture a single electron whose S = 1/2 spin interacts with the bulk excitations. 

Let us represent the impurity by a localized fermion operator Xai^)- Note that this 
fermion has no dependence upon the spatial co-ordinate r, and is a function only of the 
time T. Now we can couple this fermion to the bulk excitations which were described by C 



in Eq. (2.27) in the vicinity of the semi-metal to antiferromagnetic insulator transition. The 
full system is described by the Lagrangian C + Ci^ap where 

Amp = -hxi^pXpV^'-ir = 0,r); (2.29) 

note that whereas the Lagrangian £ is integrated over spacetime, the Lagrangian Cimp is 
only integrated over time. There are many possible couplings between the impurity and 
bulk degrees of freedom which are allowed by the symmetry of the problem, but we have 
only included a single one. This is easily seen to be the only term which is relevant under 
the RG which applies in the vicinity of the bulk quantum phase transition. 

The RG flow of the bulk-impurity coupling h was described in Refs. [2T14231 It was found 
that h approached a fixed-point coupling in the infrared, just like the couplings u and A in 
C. Thus no new couplings are necessary to describe the low energy physics of the impurity 
provided we are not too far from the semi-metal to antiferromagnetic insulator quantum 
critical point. 

Further details of the impurity dynamics may be found in the companion revieW^, where 
it was described for a closely related bulk quantum critical point. A close analogy was also 
drawn between these impurity problems and defects in super Yang-Mills theories; the latter 
can be solved^^ by intersecting brane models in string theory, and led to a description of the 
impurity criticality using a AdS2 geometry. 



F. Electrical transport 

We now revert to the honeycomb lattice without impurities. 

An important set of observables which do not acquire anomalous dimensions at the 
quantum critical point are the currents associated with global conservation laws. As the 
simplest example here, let us consider correlations of the conserved electric charge of the 
electrons, and the associated electrical conductivity a. At zero temperature (T = 0), we 
have 0" = in the insulator, while the semi-metal and the quantum critical point have finite 
non-zero values of cr, as we will now see. 

The conserved electrical current is 

= -iC^^C. (2.30) 
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Let us compute its two-point correlator, K^y{k) at a spacetime momentum k^. At leading 
order, this is given by a one fermion loop diagram which evaluates to 



d^p Tr [7^(27aPa + mp''a'')'yi,{i'ys{ks + Ps) + mp^'a'')] 



87r3 



(p2 + m'^){{p + ky + 



2 



k^ky 



dx 



x] 



\/vn? + k'^x(\ — x) ' 



(2.31) 



where the mass m = in the semi-metal and at the quantum critical point, while m = |AA^o| 
in the insulator. Note that the current correlation is purely transverse, and this follows from 
the requirement of current conservation 



h K 



0. 



(2.32) 



Of particular interest to us is the Kqq component, after analytic continuation to Minkowski 
space where the spacetime momentum fc^ is replaced by {u, k). The conductivity is obtained 
from this correlator via the Kubo formula 



a{u) = lim -^Koo{uj, k). 



(2.33) 



In the insulator, where m > 0, analysis of the integrand in Eq. (2.31) shows that that 
the spectral weight of the density correlator has a gap of 2m at A; = 0, and the conductivity 



in Eq. (2.33) vanishes. These properties are as expected in any insulator. 

In the metal, and at the critical point, where m = 0, the fermionic spectrum is gapless, and 



so is that of the charge correlator. The density correlator in Eq. (2.31) and the conductivity 



in Eq. (2.33) evaluate to the simple universal results 



Kqo{u, k) 



k' 



criu 



1 

4' 



(2.34) 



How about beyond the one-loop results? The insulator maintains a gap to charged 
excitations, and so the conductivity remains at zero. In the semi-metal, the fermions are 
gapless, but they couple only to the gapped fluctuations of the Neel order ip"-. Examination 
of the perturbation theory shows that these have no effect on the current correlators at small 



momenta and frequency, and so the results in Eq. (2.34) are exact in the limit of small uj 
and k in the semi-metal. 

At the quantum critical point, we have to consider the strong critical fluctuations asso- 
ciated with fixed-point values of the Yukawa coupling A and the quartic bosonic interaction 
u. These can be examined in the (3 — d) or the 1/A^ expansion, and require evaulation of 
multi-loop diagrams. We will not describe the computations here, but note a remarkable 
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feature: all divergences associated with the critical fluctuations cancel, and the final result 



is universal. The values of none of the couplings of the Lagrangian in Eq. (2.27) matters 



because these are all pinned by the RG fixed point. There are no anomalous dimensions, 



and the results in Eq. (2.34) generalize to 



Koo{u,k) = /C 
cr^u) = /C, 



(2.35) 



where /C is a universal number dependent only upon the universality class of the quantum 



critical point. The value of the /C for the Gross-Neveu- Yukawa model in Eq. (2.27) is not 
known exactly, but can be estimated by computations in the {3 — d) or expansions. 



III. NON-ZERO TEMPERATURES AND THE ADS/CFT CORRESPONDENCE 

We begin by some general remarks on the infiuence of a non-zero temperature, T, on a 
generic, strongly-coupled quantum critical point. Let us consider a quantum-critical point 



which has only a single relevant perturbation, s, as is the case for the model in Eq. (2.27) 
(the generalization to several relevant perturbations is immediate). So near the quantum 
critical point, the RG fiow is described by 

-r. = -s- 3.1 
at V 

In standard condensed matter notation, the eigenvalue of the relevant fiow is written in 
terms of z/, the correlation length exponent. Now in the quantum field theory in Euclidean 
time, a non-zero T corresponds to placing the theory on a cylinder of circumference 1/T. 
Such a finite size is clearly relevant in the infrared, and also indicates that 1/T transforms 
just like the temporal length under the RG. We write this as 

where z is the dynamic critical exponent. All the theories for the honeycomb lattice at half 
filling have z = 1, but we allow z to be arbitrary here. 



Eqs. (3.1) and (3.2) are of course trivial to integrate 

s{i) = se^l"" , T(£) = Te'\ (3.3) 

but the results teach us an important lesson which is summarized in the canonical quantum- 
critical phase diagram shown in Fig. |6j We ask the question: which of the relevant pertur- 
bations, s or T, is more important? To answer this question, we integrate the RG equations 
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FIG. 6: Canonical quantum critical crossover phase diagram. The dashed lines occur for T ~ l^l^*^, 
and indicate crossovers between the orange- and blue-shaded regions which are generic from any 
strongly-coupled quantum critical point. Specific features of the blue-shaded regions for the theory 



Eq. (2.27) of the transition from the semi- metal to the Neel-ordered insulator are also indicated. 
The Neel order vanishes for any T > because non-Abelian continuous symmetries cannot be 
broken in two spatial dimensions. 



to a scale (i = (!* until the winner reaches a value of unity. This allows us to conclude that 
for T > Isl^*^, thermal effects are more important than any deviation of the coupling from 
the RG fixed point. Conversely, for T < \sY'^ the couplings flow far from the critical fixed 
point before any thermally excited states need be considered. These considerations lead to 
the two distinct regimes show in Fig. [6j 

In the blue-colored regimes of Fig. [61 where T < \s\^'^ ^ the physics of the two non-critical 



phases dominates. For the model of Eq. (2.27), these are the semi-metal or antiferromagnetic 
insulator phases, both of which have well-defined quasiparticle excitations. Consequently, 
the long-time dynamics can be written using quasi-classical models of the interactions of 
these quasiparticles. 

Our interest here is primarily in the orange-colored regime of quantum criticality, T > 
\s\^'^ . Here T is the primary perturbation to the quantum critical theory. The deviation 
of the couplings from T = RG fixed point is unimportant, and the system behaves as 
it is described by the universal quantum-critical Lagrangian in the entire regime. For the 
relativistic model considered here, the strongly coupled CFT describes the dynamics of the 
orange- colored region. 

It has been arguecP^ that a central general property of quantum critical dynamics is 
the short time over which the system relaxes back to thermal equilibrium. We imagine 
perturbing the system away from equilibrium, and measuring the time, over which it 
relaxes back to local equilibrium (the adjective 'local' implies that we exclude diffusion of 
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globally conserved charges which can take a long time to reach equilibrium across the entire 
system). In the regime of strongly-coupled quantum criticality this is given by 



r„, = C^ (3,4) 

where C is a universal number dependent only upon the universality class of the transition, 
and precise definition used for Tgq. Furthermore, in all other regimes, the value of r^q is 



parametrically larger than the value in Eq. (3.4). Thus quantum criticality is described by 



a quantum fluid with the shortest possible thermal equilibration time. This characteristic 
makes it a "nearly perfect" fluid. 

It is important to note that our discussion above does not apply to CFTs in 1+1 dimen- 
sions. These are integrable systems, whose long-time dynamics is non-generic and does not 
generalize to higher dimensions. 



A. Quantum critical transport 

Let us explore the ideas above by examining the behavior of the electron conductivity 



of the model in Eq. (2.27) in the quantum-critical regime. At one-loop order, we can set 
171 = 0, and then repeat the computation in Eq. (2.31 ) at T > 0. This only requires replacing 



the integral over the loop frequency by a summation over the Matsubara frequencies, which 



are quantized by odd multiples of ttT. Such a computation, via Eq. (2.33) leads to the 
conductivity^ 



Re[(r{uj)] 
lm[(r{ijj)] 



(2Tln2)5(a;) + ^tanh('M 



°° dfi^ fRe[a{Q)] - 1/4 



TT 



u — Q 



(3.5) 



where V is the principal part. Note that this reduces to Eq. (2.34) in the limit to ^ T. 



However, the most important new feature of Eq. (3.5) arises for to <ti T, where we find a 



delta function at zero frequency in the real part. Thus the d.c. conductivity is infinite at 
this order, arising from the coUisionless transport of thermally excited carriers. 



Exercise: Evaluate i^oo from Eq. (2.31) at T > 0. First perform the trace over the 



Dirac matrices, and then the summation over the frequency. Subtract from your answer the 
result of integrating over the frequency; this subtraction can be compensated by the T = 



result in Eq. (2.31). The remaining expressions are explicitly convergent in the ultraviolet. 



and the integration over spatial momenta can be evaluated. Finally, analytically continue 



the answer to real frequencies to obtain Eq. (3.5). 
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The relaxational processes associated with Eq. (3.4) should lead to collisions between 



the thermally excited carriers and broaden the delta function at zero frequency. However, 
this relaxation does appear in a direct perturbative analysis of the critical theory in powers 
of (3 — d) or As has been discussed elsewher^^^*^, an infinite order resummation is 

required, whose simplest realization requires solution of a quantum Boltzmann equation. 
Such a solution shows that the delta function acquires a width of order (3 — d)'^T or T/N, 
and so there is a large d.c. conductivity of order (3 — c?)"^ or A^. Thus cr{u)) has the form 



of 'Drude peak' at zero frequency, and the behavior in Eq. (2.35) for to ^ T. However, the 



accuracy of such a Boltzmann equation computation is untested, and it is hkely that such 
perturbative analyses of quantum-critical dynamics are quantitatively unreliable. 

We will be satisfied here by scaling arguments which generalize the T = quantum-critical 



results in Eq. (2.35) to the T > quantum critical region in Fig. 6 The quantum-critical 



relaxational processes invalidate the form in Eq. (2.35) for the density correlation function. 



and we instead expect the form dictated by the hydrodynamic diffusion of charge. Thus for 
Kqq, Eq. (2.35) applies only for uj ^ T, while 



Koo{uj,k) = X 



Dk^ 



uj<t:T. 



lUJ 



(3.6) 



Here x is the charge susceptibility (here it is the compressibility), and D is the charge 
diffusion constant. Associated with Eq. ( |3.4 ), these have universal values in the quantum 
critical region: 

X = C^T , D = ^, (3.7) 

where again and Cd are universal numbers. For the conductivity, we expect a crossover 
from the collisionless critical dynamics at frequencies w ^ T, to a hydrodynamic collision- 
dominated form for ^ T. This entire crossover is universal, and is described by a universal 
crossover function 

lC,{u/T). (3.8) 



cr w 



The result in Eq. (2.35) applies for u ^ T, and so 



/Co-(c)o) = /C. 



(3.9) 



For the hydrodynamic transport, we apply the Kubo formula in Eq. (2.33) to Eq. (3.6) and 
obtain 

/C,(0) = C^Cd (3.10) 



which is a version of Einstein's relation for Brownian motion. 
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B. The AdS/CFT correspondence 



Portions of this section have been adapted from Chapter 15 of Ref. [191 
It turns out the AdS/CFT correspondence is 'just what the doctor ordered' to compute 
strongly-coupled quantum critical dynamics and transport in the orange-colored region of 
Fig. [6j This is a consequence of a crucial property: even at the level of the classical grav- 
ity approximation in the AdS theory, the system relaxes back to thermal equilibrium in 



a time which obeys Eq. (3.4). No other method in condensed matter physics shares this 
remarkable feature. We will review specific computations by this method of the univer- 
sal function lC„{uj/T), and of the coUisionless-to -hydrodynamic crossover in the density 
correlation function. 

The CFT solvable by the AdS/CFT correspondence may be viewed as a generalization of 



the CFT described by Eq. (2.27). It has a closer resemblance to the SU(2) gauge theory we 
consider later in Eq. (5.8). We take the structure of critical matter fields coupled to a gauge 
field, and generalize it to a relativistically invariant model with a non-Abelian SU(A^) gauge 
group and the maximal possible supersymmetry. The resulting supersymmetric Yang-Mills 
(SYM) theory has only one independent coupling constant g. Under the RG, it is believed 
that g flows to an attractive fixed point at a non-zero coupling g = g*', the fixed point 
then defines a supersymmetric conformal field theory in 2-1-1 dimensions (a SCFT3). We 
are interested here in computing the transport properties of the SCFT, as a paradigm of 
quantum critical transport at a strongly interacting quantum critical point. 

The solution proceeds by a dual formulation as a four-dimensional supergravity theory 
on a spacetime with uniform negative curvature: anti-de Sitter space, or AdS4. The solu- 
tion is also easily extended to non-zero temperatures, and allows direct computation of the 
correlators of conserved charges in real time. At T > a black hole appears in the gravity 
theory, resulting in an AdS-Schwarzschild spacetime, and T is the Hawking temperature of 
the black hole; the real time solutions also extend to T > 0. 

The reader is referred to the original papeil^, and to the TASI lectures by Son for an 
explicit description of the method. In the AdS/CFT correspondence, every globally con- 
served quantity in the CFT gets mapped onto a gauge field in AdS. Moreover, in the leading 
classical gravity theory on AdS, different global charges commute with each other, and so 
can be considered separately. In the end, we have a U(l) gauge field on AdS for every global 
conservation law of the CFT. The low energy effective field theory on AdS4 has the standard 
Einstein-Maxwell action for gravity-|-electromagnetism: 



94 



4 



(3.11) 



Here gab is the AdS-Schwarzschild metric {g is its determinant). Fab is the Maxwell flux 
tensor, and g4 is a dimensionless coupling constant fixed by the value of N in the SU(A^) SYM 
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FIG. 7: Spectral weight of the density correlation function of the SCFT3 with M = 8 supersym- 
metry in the collisionless regime 
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FIG. 8: As in Fig. [7| but for the collision- dominated regime. 



theory. This 4-dimensional Maxwell theory can be used to compute the density correlation 
function, Koo(uj,k), of the 3-dimensional SYM theory, and the results are shown in Fig. [T] 
and |8| The most important feature of these results is that the expected limiting forms 
in the collisionless (Eq. (2.35)) and collision-dominated (Eq. (3.6)) are obeyed. Thus the 
results do display the collisionless to collision-dominated crossover at a frequency of order 



ksT/h, as we expected from the physical discussion in Section III A 



At this point, we describe some technical aspects of the results which turn out to have 
important physical implications. For this, let us generalize the constraints on K^^, from 
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current conservation in Eq. (2.32) to non-zero temperatures. At T > 0, we do not expect 



K^i, to be relativistically covariant, and so can only constrain it by spatial isotropy and 
density conservation. These two constraints, along with dimensional analyses, lead to the 
most general form 



K^,{uj, k) = VW^{P^^ K^iu, k) + P^, K\uj, k) 



(3.12) 



where = {—uj,k) and k = \k\. The K^'^ are dimensionless functions of the arguments, 
and depend upon u and the magnitude of the 2- vector k. Also and Pf^, are orthogonal 



projectors defined by 



P'^ 

-'00 



P^ 







ki kj 



k^ky 



P. 



(3.13) 



with rjfj^y = diag(— 1, 1, 1), and the indices i,j running over the 2 spatial components. Thus, 
in the general case at T > 0, the full density and current responses are described in terms 
of two functions K^''^{k, u), representing current fluctuations longitudinal and transverse to 
the momentum. These two functions are not entirely independent. At T > 0, we expect all 
correlations to be smooth functions at /c = 0: this is because all correlations are expected 
to decay exponentially to zero as a function of spatial separation. However, this is only 



possible from (3.12) if we have the additional relation 



(3.14) 



Finally, we note that application of the Kubo formula in Eq. (2.33) to Eq. (3.12) yields 



(3.15) 



The relations of the previous paragraph are completely general and apply to any theory. 
Specializing to the AdS-Schwarzschild solution of SYM3 as determined by the Einstein- 



Maxwell theory in Eq. (3.11), the results were found to obey a simple and remarkable 
identit}^ 



K^{u,k)K^{u,k) =lC 



(3.16) 



where is a known pure number, independent of u and k. This identity is a consequence of 



the self-dual structure of Eq. (3.11): the Maxwell action is in 3+1 dimensions is well-known 



to have a self-dual structure corresponding to the exchange of electric and magnetic fields. 
Thus we have the important and key result that every global charge in a CFT3 maps onto 



a self-dual theory in the leading gravity approximation on AdS4. The identity in (3.16) is a 
consequence of this emergent self-duality of CFT3s. 



The combination of (3.16) and (3.14) now fully determine the response functions at zero 
momenta: K^{u), 0) = K'^{uj, 0) = /C. Computing the conductivity from Eq. (3.15), we then 
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have 



cr(a;) = }Ca{ou/T) = /C; 



(3.17) 



I.e. 



the scahng function in Eq. (3.8) is independent of co and equal to the value in Eq. (3.9). 



This result is an important surprise: the conductivity of the classical gravity theory on AdS4 
is frequency-independent. Furthermore, its value is fixed by self-duality to be the constant 



/C appearing in the self-duality relation (3.16). All these remarkable results are a direct 



consequence of the self-duality of the U(l) Maxwell theory on AdS4. 



Given the strong consequences of self-duality relation in Eq. (3.16), it is useful to ask 



whether it can be valid for CFTs beyond those described by the classical Einstein-Maxwell 
theory on AdS4. This question was addressed recently by Myers et alJ^. The examined the 



general structure of the higher-derivative corrections to Eq. (3.11), and argued that for the 



current correlations the leading terms could always be transformed into the following form 
which has only one dimensionless constant 7 [L is the radius of AdS4): 



M 



1 

7, 



d Xy 



-9 



ab j-^cd 



(3.18) 



where the extra four-derivative interaction is expressed in terms of the Weyl tensor Cabcd- 
A crucual observation of Ref . 1551 was that stability and causality constraints on the effective 



theory restrict I7I < 1/12. A generalized duahty relation applies also to S'^^j. However this is 
not a se//-duality. The dual GET has current correlation functions which were characterized 
by functions K^'^{u, k) which are distinct from those of the direct GET K^'^{uj, k), and the 



self-duality relation of Eq. (3.16) take the less restrictive form 



K^{uj,k)K^{Lo,k) = K? 



K^iuj,k)K^{uj,k) = K? 



(3.19) 



These duality relation determines the correlators of the dual GET in terms of the direct 
GET, but do not fix the latter. Determination of the functions K^{u, 0) = K'^[u, 0) requires 
explicit computation using the extended theory Syec, and the results for the conductivity are 
presented in Eig. |9] Now the conductivity does have a non-trivial universal dependence on 
Lo/T. However, as is clear from Eig, [9| stability conditions on the effective theory on AdS4 
allow only a limited range of dependence on tu/T. The smooth tu/T dependence in Eig. [9] 
should be contrasted to the very singular dependence in the free-field result in Eq. ( |3.5[ ); 
the former is clearly more generic for a strongly-coupled GET. It is also interesting to note 
that the u/T dependence in Eig. [9] for 7 > is very similar to the structure we discussed 



in Section HI A on the effect of collisional broadening of the singularities in Eq. (3.5): the 



AdS4 result shows a collision-dominated Drude peak at w = 0, and a coUisionless critical 
continuum at large u. 

Does the duality mapping of Myers et al^ have an interpretation directly in the GET, 
without using the mapping to AdS4? It has been arguecP^^ that this duality is the analog 
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FIG. 9: Frequency dependent conductivity for CFT3s for which the AdS4 theory includes the lead- 
ing correction beyond the Einstein-Maxwell theory from Ref. I33l The co-efficient of this correction 
in the action is 7, and stability requirements impose the bound I7I < 1/12. 



of the 'particle-vortex' duality of condensed matter physics. The latter is an exact duality 
of the critical theory of a complex relativistic scalar fielcP^ (this the theory in Eq. (2.27) 



without the fermions, and with ip"' having two components). In the particle- vortex duality, 
the world line of the complex scalars is reinterpreted as the world line of vortices in the 
dual theory of a dual complex scalar interacting with an emergent electromagnetic field. 
This particle-vortex duality also allows us to interpret the structure of the results in Fig. [9] 



for 7 < 0. Note that from Eq. (3.19) the conductivity of the direct CFT maps onto the 
resistivity of the dual CFT. Thus for 7 < 0, it is the dual CFT which will have a conductivity 
which has the structure discussed in Section [TlIA[ with a collisionally-broadened Drude peak 



at w = 0. Thus we conclude that a Boltzmann-like picture of transport applies better to the 
particle interpretation of the CFT for 7 > 0, and to the vortex interpretation for 7 < 0. 

Let us summarize the lessons we have learnt from the AdS theory of quantum critical 
transport in strongly interacting systems in 2+1 dimensions. This theory should be view 
as complementary to the quasiparticle-based theory, whose implications were discussed in 
Section IIII A[ The lessons are: 

• There are a large class of strongly interacting 2+1 dimensional quantum liquids which 
relax back to thermal equilibrium in the shortest possible time of order h/{kBT), as 



we indicated in (3.4). 



The quasiparticle transport theorjJ^SlE] starts from the free theory with an infinite 
thermal equilibration time, and includes the effect of weak interactions using the 
Boltzmann equation. Complementary to this is the quantum-critical transport theory 
applicable for the shortest possible equilibration time of order h/{kBT), which is the 
classical Einstein-Maxwell theory on AdS4. 
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• The Einstein-Maxwell theory exhibits collisionless dynamics for u ^ T, and collision- 
dominated dynamics for u <t^T, as we displayed in Figs. [7] and |8j 

• All continuous global symmetries are represented by a self-dual Einstein-Maxwell the- 
ory. 

• This emergent self-duality implies that, in systems with particle-hole symmetry, (j{uj) 
is frequency-independent in the Einstein-Maxwell theory and equal to the self-dual 
value. 

• For systems with particle-hole symmetry, a frequency dependent conductivity is 
obtained^ upon considering corrections to the effective Einstein-Maxwell theory, with 
the forms in Fig. |9j Stability conditions on the effective theory allow only a limited 
range of frequency dependence, which depend upon the single parameter I7I < 1/12. 
For 7 > 0, the frequency dependence has the form expected from collision-dominated 
transport of particles, while for 7 < it is characteristic of the transport of vortices. 
It is remarkable that the physical pictures expected from the Boltzmann transport 
analysis correspond to precisely those obtained from the holographic approach. 

• Such quantum-critical fluids also have universal momentum transport. By extending 
the scaling arguments to momentum transport we would conclude that the ratio of 
the shear viscosity to the entropy density rj/s should equal a universal number charac- 
terizing the collision-dominated regime. This number was computed in the Einstein- 
Maxwell theory by Kovtun et a/.'2SE3 and found to equal h/^ATiks)- The shortest 
possible relaxation time implies that 77 is also the smallest possible and so these fluids 
are 'nearly perfect'. 



IV. U(l) GAUGE THEORY AND THE VALENCE BOND SOLID ON THE HON- 
EYCOMB LATTICE 



We now return to the honeycomb lattice at half-filling. In Section II D we described a 
quantum phase transition in which two characteristics of the ground state changed simul- 
taneously. In the charge sector, the one electron excitation gap opened up leading to a 
transition from the semi-metal to the insulator. And in the spin sector, the breaking of 
SU(2) spin rotation symmetry led to Neel order in the insulator. However, in the Mott 
picture, the insulating behavior is tied to repulsion between the electrons, which keeps them 
apart, rather than to any specific symmetry breaking. This would suggest that it is possible 
to have an insulating state while preserving spin rotation invariance. We will explore such 
a possibility in the present section. 

Readers not interested in issues related to electron fractionalization and emergent gauge 



fields in insulators may skip ahead to the discussion of metallic phases in Section VI 
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Our approach will be begin in the Neel-ordered insulator, and restore spin rotation in- 
variance by allowing for slow angular fluctuations in the local orientation of the Neel order 
parameter At the same time, it will also pay to transform the fermions to a 'rotating 
reference frame' so that their spin is measured relative to the orientation of the local Neel 
ordei'^Smni, This transformation is most conveniently done using spinor variables. So let us 
decompose the vector Neel order (yj" into a complex two-component bosonic spinor Za by 

^'' = z:<0zp (4.1) 



Such a decomposition is familiar from early work of D'Adda et al!^ and Witten!^, who 
established the equivalence between the 0(3) non-linear a— model and the CP^ model in 2 
spacetime dimensions. A similar equivalence does not immediately apply in the 3 spacetime 
dimensional case of interest here because point defects in spacetime have to be treated with 
some card^. In particular, the theory for the fluctuations of the vector field must allow for 
point spacetime defects ('instantons') where ip°' = 0, which are known in the condensed mat- 
ter literature as 'hedgehogs'. Note that these hedgehogs are present even in a 'fixed-length', 
non-linear cr-model in which we set ip"''^ = 1; such models require ultraviolet regularization, 
and the hedgehogs are invariably permitted in the regulated theory e.g. with a lattice reg- 
ularization. Ignoring these defects momentarily, let us proceed as in the earlier workP^El_ 



The parameterization in Eq. (4.1) is invariant under the U(l) gauge transformation 

Za ZaC'^ (4.2) 

and so the theory for the Za must be a U(l) gauge theory involving an emergent U(l) gauge 
field Af^. The boson only terms in Eq. (2.27) are equivalenlP^'^ to a U(l) gauge theory for 



the complex scalars Za 

= \{d^- iA^)z^\^ + s\z^\^ + u{\z^\^f. (4.3) 

Here the gauge field is dynamical, and will acquire a Maxwell action after high energy 
Za modes have been integrated out. 



Let us now discuss the point defects. Eqn (4.1 ) implies^that the hedgehogs in ip^ become 



Dirac monopoles in A^: these are tunnelling events associated with a change in the total A^ 



flux by 27r. Such monopoles are permitted by the U(l) gauge theory in Eq. (4.3) only if the 
U(l) gauge field is compact. So we must account for the dynamics of such a compact U(l) 
gauge theory to completely account for the fluctuations of the local antiferromagnetic order. 
The dynamics of the matter fields can suppress the monopoles in some case^^^, and this 
can then lead to deconfined critical points or phases with a gapless U(l) photon excitation 
associated with an effectively non-compact U(l) gauge field. We will find an example of this 



phenomenon in Section VII 
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Let us now turn to the fermionic excitations in antiferromagnetic insulator. We transform 
these to the rotating reference frame by writin^^BHlQ] 



(4.4) 





where ipp^ p = ±, are the "electrons" in the rotating reference frame. The index p measures 
the spin-projection along the direction of the local Neel order. However, more properly it 



is the "charge" under the emergent U(l) gauge field because Eq. (4.4) is invariant under 

(4.5) 



Eq. (4.2) and 



An important consequence of Eqs. (4.4) and (4.1) is the identity 



\2\2 



(4.6) 



Thus the effective moment acting on the ijj fermions is always along the expected 
by the transformation to a rotating reference frame. 

Let us now take the continuum limit for the fermions in the rotating reference frame. We 



follow exactly the same mapping as in Eq. (2.8) to map the lattice if) fermions to continuum 



8-component ^ fermions. Based upon Eq. (4.6), we also expect the ^ fermions to experience 



a field polarized along the direction. Combined with gauge invariance, and the structure 



of Eq. (2.27), we are led to the following Lagrangian density for 



(4.7) 



Note this is the theory of Dirac fermions of mass \XNq\ coupled to a U(l) gauge field. The 



coupling of the fermions to the U(l) gauge field in Eq. (4.7) can also be derived explicitly by 



substituting Eq. (4.4) into the last term in Eq. (2.27), and using the expression for the U(l) 



gauge field in the CP^ model. We will present a more explicit derivation of an emergent 



gauge field for the case of the triangular lattice in Section |Vl] below. 

The Lagrangian + is then our U(l) gauge theory for the fluctuating Neel state. 



complementary to the Gross-Neveu- Yukawa theory in Eq. (2.27). The remainder of this 



section is devoted to understanding the physical properties of + d^i. 
Let us now discuss the phases of this U(l) gauge theory. 

First, we have the Higgs phase, where s < and is condensed. Here the U(l) photon 
is gapped, and spin rotation invariance is broken. This is just the insulating Neel state, and 



its properties are identical to the Neel ordered state described by Eq. (2.27). 



The other phase with s > has Za gapped and spin rotation invariance is preserved. 



However, as is clear from Eq. (4.7), the fermionic spectrum remains gapped. Thus this 



phase is clearly not the semi-metal of Eq. (2.27). Instead it is a new insulating phase with 
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spin rotation invariance preserved. Thus we have achieved our objective of describing an 
insulator without Neel order. 

However, this insulator is not a featureless state with a spin and a charge gap, as we will 
now show. The interesting physics arises from an interplay of the monopole events with the 
filled band of fermionic states. If we integrate out this filled band via £4,, we generate an 
effective Maxwell action for the U(l) gauge field 

We recall that the U(l) gauge field is compact, and it was shown by Polyako'v'^ that such 
a gauge field always acquires a mass gap and is in a confining phase in 2+1 dimensions. 
The confinement is caused by the proliferaction of monopole tunneling events. Here we will 
sho'ftff^E^ESEH] lY^Q^ f^]^Q monopole operator has non-trivial transformation properties under 
the symmetry group of the honeycomb lattice: consequently, the proliferation of monopoles 
in the confining phase breaks a lattice symmetry due to the appearance of valence bond 
solid (VBS) order. 



To see this, it is useful to add an external source to the fermion Lagrangian in Eq. (4.7) 
so that £^ becomes 

U = ^lf.{d^ + tcr'A^)^ - AiVo^pV^^ - '-B^^^^p'^. (4.9) 

This source has been judiciously chosen so that when we integrate out the fermions, the 
action for the gauge field in Eq. (4.8) acquires a mutual Chern-Simons termpEHl 



= 10 u (W^'--^^)' + ^B^e^^xd.A^ (4.10) 

i/TT AiVn 27r 



Let us now proceed with Polyakov's duality mapping on Eq. (4.10) to obtain an effective 
theory of monopoles: the Bfj_ source term will allow us to deduce the connection between 
the monopole operator and the underlying lattice fermions. The first step corresponds to 
decoupling the Maxwell term by a Hubbard- Stratonovich field, Y^, to obtain 

'^A = + ^Yf^^f^uxduAx + +^B^e^,xduAx (4.11) 

Now, we integrate over A^, and this yields the constraint 

Y^ = d^(l)-B^. (4.12) 

where (p is the scalar field which is dual to the photon. We have judiciously chosen factors of 
(27r) above to ensure a normalization so that e*"^ is the monopole operator. Finally, inserting 
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Eq. (|4.12|) into (|4.11[) we obtaiiPI^Hl 



3|AAr, 
47r 



o| 



(4.13) 



In the absence of the external source this is a free scalar field theory, which implies that 
the monopole operator e*"^ has long-range correlations in 2+1 dimensions. In other words, 



the free photon phase described by Eq. (4.8) has a non-zero VEV with (e*"^) 7^ 0. 



The term in Eq. (4.13) will help us link the monopole operator to the underlying 
electron^. First, we notice that the theory in Eq. (4.9) actually enjoys a gauge invariance 
under which 



\E' — )• exp 



4. 



Bp, — )• Bp 



dp9 



(4.14) 



where ^ is a field with an arbitrary spacetime dependence. (Note that this gauge invariance 



is distinct from that associated with the gauge field in Eq. (4.5), under which \E' — )■ 
exp(— 'icr^(^)\E'.) Now we observe that this gauge invariance extends also to Eq. (4.13), under 
which 

^ e^^e*"^. (4.15) 



The combination of Eqs. (|4.14l) and (|4.15l) now allows us to identify the operator e*"^. We 



look for a fermion bilinear of the form \E'M\I' so it transforms like Eq. (4.15) under the gauge 



transformation in Eq. (4.14). Moreover, the Lorentz invariance of the theory implies that 



the matrix M should commute with the 7^ matrices in Eq. (2.12). This leads us to the 
unique choice 

e'^ ~ Wr^(p^ + ipy)m ~ Cr^(p^ + ipy)C (4.16) 



It now remains to use the geometric definitions in Section II A and Eq. (2.8) to deduce the 



physical interpretation of the fermion bilinear in Eq. (4.16). A careful analysi^ along these 



lines shows that e**^ is an operator associated with the valence bond solid (VBS) order in 



Fig. 10, and the VEVs of the operators in Eq. (4.16) imply long-range VBS order. 



Exercise: Compute the transformations of the fermion bilinear in Eq. (4.16) under 



honeycomb lattice symmetries such as translations, refiections, and rotations by 60 degrees. 



All these transformations map the pattern in Fig. 10 either to itself or to 2 equivalents 
patterns. Assign the weights 1, e^'^*/^, and e^'^*/'^ to these patterns, and show that their 



transformations coincide with those of Eq. (4.16). 



Thus we reach our main conclusion: the insulating phase without Neel order as described 
by the U(l) gauge theory Cz + C-^, has long-range VBS order. This order onsets with the 
confinement induced by the proliferation of monopoles. 

It is interesting to note that the matrices in the fermion bilinears associated with VBS 
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FIG. 10: A schematic illustration of the valence bond solid (VBS). The ellipses represent singlet 
valence bonds between the electrons. These reside preferentially in the pattern shown in the VBS 
state. Expectation values of all spin-singlet observables, such as SfSj or c\^Cja + c^j^^^Cia, are 
different on the links with the ellipses than those without. 



order ~ Ct'^{p^ + ip'^)C and Neel order ~ Cp^a"'C all anti-commute with each other, and 
commute with the 7^ matrices in Eq. (2.12). This can be used to formulate the present 
theory without gauge fields but using a Wess-Zumino-Witten termpHSZl -^q -^[W ^ot explore 
this approach here. 

We can also use the methods of this section to address the nature of the transition between 
the Neel and VBS phases. We will not go into details here, but this transition has been 
proposecP^'^ to be a deconfined critical point at which monopoles are suppressed, and the 
critical theory is the non-compact version of the U(l) gauge theory given by C^- 



V. SU(2) GAUGE THEORY AND A PHASE DIAGRAM FOR THE HALF- 
FILLED HONEYCOMB LATTICE 



Sections [IT] and IV have so far described 3 possible phases of the honeycomb lattice at 
half-filling: the semi-metal, the insulator with Neel order, and the insulator with VBS order. 
The first two phases appear in the theory C in Eq. (2.27), while the latter two appear in the 
theory in Eqs. (4.3) and (4.7). This is an unsatisfactory state of affairs: we would like 

to write down a single unified theory in which all 3 phases appear. One approach, implicitly 
mentioned above, is to extend Eq. (2.27) by including an additional two-component real 
scalar field representing the VBS order parameter, and couple it to the fermion bilinear 
Cr^(p^ + ip'^)C appearing in Eq. (4.16). Integrating out the fermions in the background 
of a spatially varying 5-component scalar representing the Neel and VBS orders yields a 
Wess-Zumino-Witten term for the scalar fielcl^SMS] rpj^g resulting theory is difficult to work 
with, and little is known about it in the regime where all three phases can meet. 

Here we will present an alternative approach which allows for exotic phases using an 
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emergent SU(2) gauge field. We will find that the resulting phase diagram has a fourth 
semi-metallic phase with an emergent topological order, and an interesting multicritical 
point. 

Our starting point is the observation that the decomposition of the electron into spinful 



bosons and spinless charged fermions in Eq. (4.4) has a larger gauge invariance^ than U(l) 



Rewriting eq. (4.4) using a natural matrix notation 



c = 



where 



R = 



(5.1) 



(5.2) 



we note that Eq. (5.1) is invariant under the gauge transformation generated by SU(2)g 



matrix U under which 



R 



c ^ c. 



(5.3) 



This SU(2)g gauge transformation should be distinguished from the global SU(2) spin rota- 
tion V , under which 

R^VR , , c^Vc. (5.4) 

Turning to the Neel order y?", this clearly transforms as a 3 under the global SU(2). 



However, the parameterization for the Neel order in Eq. (4.1) is not invariant the SU(2) 



gauge tranformation in Eq. (5.3). As written, Eq. (4.1) is invariant only under the U(l) 



gauge transformation in Eq. (4.2) which was the reason for our original choice of a U(l) 



gauge theory in Section IV Thus we cannot use Eq. (4.1) as our definition of the Neel order 



in the present SU(2) gauge theory. 

It is more natural to proceed herd^^EG! i^y defining the scalar fields using bilinears of the 



fermions. Thus using Eq. (2.25) and extending to continuum 8-component fermions near 
the Fermi points, we define 

ip^ = Cp'a^C. (5.5) 

From this definition it is clear that (p"" transforms as a 3 under the global SU(2), while it is 
invariant under SU(2)g, just as expected. 

We can also define the corresponding scalar using the ijj fermion^^ 



(5.6) 



Now transforms as a 3 under the gauge SU(2)g, while it is invariant under SU(2). 
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From Eqs. (5.1), (5.5) and (5.6), we find that the scalar fields are related by 



(5.7) 



These relations generalize Eq. (4.1) from the U(l) gauge theory. 

To summarize, the matter fields of our SU(2)g gauge theory are the bosonic matrix i?. 



the fermions ^, and the scalar As in Section IV, we will also need an emergent dynamic 
SU(2)g gauge field A". Using symmetry and gauge invariance, we can now write down the 
following Lagrangian densitj^ for the SU^ gauge theory; this combines and generalizes the 



Gross-Neveu- Yukawa model in Eq. (2.27), and the U(l) gauge theory in Eqs. (4.3) and (4.7) 



+ 



(9^$" - 2eAfec4$^) +^'^ 
+ Tr [{d^R - iAlRa''){d^R^ + lA^a'^R^)] 
+ s Tr [RR^) + u [Tr [RR^)]" . 



a2\2 



+ i (*") 



(5.8) 



This Lagrangian combines all three phases discussed so far, and forms the basis of our 
remaining discussion of the honeycomb lattice at half-filling. 

In mean-field theory, the model Cg has 4 phases, depending upon whether one or both 
of the scalar fields and R are condensed. These 4 phases can be identified using the 



methods developed in Section II D and IV, and lead to the phase diagram shown in Fig. 11 



First, we describe how Cg reproduces the phases and phase transitions already discussed: 

• The Higgs phase where {R) ^ breaks SU(2)g completely. Using SU(2)g gauge in- 
variance we may as well set R = 1. Then from Eq. (5.7), we have ~ (p"", the Neel 



order parameter. Also the gauge boson A't is gapped and can be neglected. Then 



the theory Cg reduces to the Gross-Neveu- Yukawa model in Eq. (2.27). As discussed 



in Section II D, this theory has semi-metal and insulating Neel phases, and these are 



shown in Fig. 11 



The Higgs phase where (A^") ^ breaks SU(2)g down to U(l). Then only the (say) 
gauge boson is active, and the theory Cg reduces to the U(l) gauge theory + jC* 



discussed in Section IV The insulating Neel and insulating VBS phases found there 



are also shown in Fig. 11 



The possible new phase of Cg is the deconfined phase where both and R are gapped. 
Then the low energy theory of Cg is simply massless QCD with the Lagrangian density 



CQCD = '^l,{d, + ia''Al)^>. 



(5.9) 
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FIG. 11: Schematic phase diagram of the SU(2)j gauge theory Cg in Eq. (5.4). The two phases in 



the bottom are described by the Gross-Neveu- Yukawa model in Eq. (2.27), while the two phases 



on the left are described by the U(l) gauge theory Cz + ^C* in Eqs. (4.3) and (4.7). 



This QCD with a SU(2)g gauge group with massless 2-component Dirac fermions which 
carry Nc = 2 colors and Nf = 2 flavors. When Nf is large enough, it can be shown from a 
1/Nf expansion that the confining tendencies of non-Abelian gauge fields are screened, and 
Cqcd describes a non-trivial CFT, with anomalous dimensions for all observables which are 
not currents of global flavor or spacetime symmetries. It is an open question whether such 
a critical phase is allowed for Nf = 2, as we have assumed in Fig. 11 If not, then this 



phase will be unstable to confinement into one of the other phases of Fig. 11 If present, 
this deconfined phase would be a topologically ordered semi-metal with a spin gap; it is an 
'algebraic charge liquid' (ACL) in the notation of Ref. |65l The gapless ^ fermions carry 
electromagnetic charge, and so there is no gap to charged excitations excitations and this 
phase is not an insulator. However, the \1/ fermions are spinless, and SU(2) spin is only 
carried by the gapped bosonic excitations; hence the spin gap. 
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Fig 11^ also shows an interesting multi-critical point M, where all 4 phases meet; if the 
massless QCD phase is confining, this would be the meeting point of 3 phases. Here the 
SU(2)p gauge bosons, the scalars and i?, and the fermions \1/ are all gapless and critical. 
Thus M realizes a non-trivial CFT which can be perturbed by the two relevant directions 
of the plane of Fig. H Indeed, it is not unreasonable to view this multicritical M theory 
as a non-supersymmetric analog of the M-theory of strings! 



VI. THE METAL-INSULATOR TRANSITION ON THE TRIANGULAR LAT- 
TICE 



This section will describe possible phases of the Hubbard model in Eq. ( 1.1 ) on the trian- 
gular lattice. We will now consider the case of generic density, so that unlike Sections [H, IV 
and [V] on the honeycomb lattice we will allow (njf),(^jt) 7^ 1/2, although the half-filled 
density will also appear in our phase diagram. Unlike the honeycomb lattice, we will ig- 
nore the possibilities of magnetically ordered phases in which the global SU(2) spin rotation 
symmetry is broken. The half-filled model on the triangular lattice likely does have antifer- 
romagnetic order in the limit of large f/, but we will not consider this complication here. 
Our purpose here is to describe the structure of possible phases without magnetic order. 

The most significant difference from the honeycomb lattice is apparent in the limit of 
small f/, when the electrons are nearly free. Then the triangular lattice ground state is 
a metal at all densities, unlike the semi-metal state found on the honeycomb lattice at 
half-filling. The semi-metal had a spectrum with a relativistic structure at low energies, a 
fact which we have exploited in our discussion so far. However the metal has zero energy 
excitations along a line in momentum space, the Fermi surface, and the fermionic excitations 
near the Fermi surface do not have a relativistic spectrum. 

Landau's Fermi liquid (FL) theory provides a complete description of the universal prop- 
erties of the metal. We will not review this theory here: the reader is referred to Chapter 
18 of Ref. [19] for the author's perspective. A discussion in the context of the gauge-gravity 
duality appears in a recent papei'^. For our purposes here, we need only two basic facts: 
{i) the fermionic excitations near the Fermi surface are essentially non-interacting electrons, 
and {a) the area enclosed by the Fermi surface is equal to the electron density — this is 
Luttinger's theorem, which we state more explicitly below. 



The FL metal can be described by ignoring the U interactions, and transforming Eq. ( 1.1 ) 



to momentum space. Unlike the honeycomb lattice, there is only one site per unit cell of 



the triangular lattice, and so the analog of Eq. (2.6) is now simply 



i/o = ^cJ^(fc) — /i — 2i(:^cos(fc ■ ei) + cos(fc • 62) + cos(fc ■ 63) j Ca{k), (6.1) 



where there are no Pauli matrices associated with sublattice index, the are as in Eq. (2.1), 
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FIG. 12: The triangular lattice 




FIG. 13: The electronic dispersion in Eq. (6.1) for n = and t = 1 



and are shown in Fig. 12, The reciprocal lattice now consists of the vectors ^^UiGi, where 



Eq. (2.2) is replaced by 



Air 



ei - 62 



Go 



An 



162 - egj 



G, 



[63 - ei, 



(6.2) 



while the first Brillouin zone is a hexagon with vertices given by Eq. (2.3), as shown in 



Fig. 14 The electronic dispersion in Eq. (6.1) is plotted in Fig. 13 it only has simple 



parabolic minima at fc = 0, and its periodic images at fc = G, and there are no Dirac 
points. At any chemical potential, the negative energy states are occupied, leading to a 



Fermi surface bounding the set of occupied states, as shown in Fig. [T4j Luttinger theorem 
states that the total area of the occupied states, the shaded region of the first Brillouin zone 
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FIG. 14: The Fermi surface of Eq. (6.1) for /i = 1/2 and t = 1; the occupied states are shaded. 



Also shown are the periodic images of the Fermi surface in their respective Brillouin zones. 



in Fig. 14 occupies an area, A, given by 



27r2 



(6.3) 



where Af = J2a '^a^^a is the total electron density. This relationship is obviously true for 
free electrons simply by counting occupied states, but it also remains true for interacting 
electrons, as reviewed recently in Ref. |5l 

Now we turn up the strength of the interactions, U. For the honeycomb lattice, we pre- 



sented in Sections II C and II D a treatment which allowed for the appearance of spontaneous 
magnetic moment on each site, leading to the onset of antiferromagnetic order at large U. 
We also found that the onset of antiferromagnetic order co-incided with appearance of the 
insulator i.e. the semi-metal to insulator transition. This co-incidence was related to the 
appearance of a gap in the spectrum even for an infinitesimal antiferromagnetic moment. 



as in Eq. (2.21). We can apply a similar treatment here to the triangular lattice. However, 



such an analysis finds that the onset of magnetic order does not co-incide with the metal- 
insulator transition. Instead, we find an intermediate metallic phase with magnetic order, in 
which the Fermi surface has been reconstructed into small "pockets" . Such a reconstruction 
is observed in many correlated electron systems. However, we will not explore this route to 
the insulator here, and refer the reader to recent paperd^SEll 

Instead, we will explore here a distinct route to the destruction of the Fermi liquid, one 
which reaches directly to an insulator which is a 'spin 

liquid'SSMZI. The spin hquid insulator 
is a phase in which the spin rotation symmetry is preserved, and there is a gap to all charged 
excitations. In these respects, the spin liquid is similar to the insulating state discussed in 



Section IV However, in Section IV we found that such a insulator had an emergent U(l) 
gauge field, A^, and the proliferation of monopole defects in led to a confining phase in 
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which lattice translational symmetry was broken by the appearance of VBS order. Here we 
will find that the triangular lattice spin hquid also has an emergent U(l) gauge field, but the 
presence of Fermi surfaces in the spinful excitations leads to the suppression of monopole 
events. Consequently, we have a deconfined phase with gapless gauge excitations, no lattice 
translational symmetry breaking, and the spin liquid character survives. 

The key to the description of the metal insulator transition is exact rewriting of the 



Hubbard model in Eq. (1.1) as a compact U(l) lattice gauge theory. We proceed with a 



method which parallels that in Section IV, of transforming to a 'rotating reference frame'. 
However, instead of using the frame of reference of local antiferromagnetic order, we use a 
quantum rotor to keep track of the charge on each lattice site. Each rotor has a periodic 
angular co-ordinate di with period 27r; hence the states of the rotors are e^^"'^^ where n^j is 
a rotor angular momentum, whose eigenvalues take all positive and negative integer values. 
We will use the state with all n^i = to represent the states with one electron each lattice 



site. The analog of the transformation to a rotating reference frame in Eq. (|4.4|) is no"^ 

J a 



(6.4) 



where we have dropped the implicit site index, and fa are neutral fermions ('spinous') which 
keep track of the orientation of the electron. We can now identify the 4 states on each lattice 



site in Eq. (2.13) with corresponding states of the rotor and spinous: 



|0) 




e-^'^lO) 


4|0) 




/I|0) 


4410) 







Note that these allowed states obey the constraint 

fafa 1- 



(6.5) 



(6.6) 



Associated with this constraint is the U(l) gauge invariance which is the analog of Eqs. (4.2) 

fa^fae'"^ , ^^^ + C. (6.7) 



and (4.5) 



Just as in Section IV, there will be an emergent gauge field An in the effective theory of this 



model. The constraints in Eq. (6.6) will be the Gauss law of this theory. 

First, let us rewrite the Hubbard model in terms of these new variables. Our degrees 
of freedom are the Fermi operators fia on each lattice site which obey the usual canonical 



fermion anti-commutation relations as in Eq. (1.3), and the rotor angle 'di and angular 
momentum Uri which obey 

['&i,nrj] = i6ij. (6.8) 
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The Hubbard Hamiltonian in Eq. (1.1) is now exactly equivalent to 



H[f,^] = -J2t^jfLfjo.e 



-/i(nri + 1) + ^ririiriri + 1) 



(6.9) 



provided our attention is restricted to the set of states which obey the constraint in Eq. (6.6) 



on every lattice site; note that the Hamiltonian in Eq. (6.9) commutes with constraints in 



(6.6), and so these can be consistently imposed. In Eq. (6.9) we have used the rotor angular 



momentum to measure the charge on each site, and so the dependence of the energy on (i 
and U can be expressed in terms of rin alone. 

We can now implement the commutation relations, the Hamiltonian, and the constraint 
in a coherent state path integral 



Z 



VfUr)Vfl{r)VUr)Vn„{r)V\{T) exp | 



drH[f,^] 



dr 



V 



t dfu 



dr 



(6.10) 



where '&i{T) takes values on a circle with unit radius, ensuring quantization of eigenvalues 



of the angular momentum n^j to integer values. The constraint in Eq. (6.6) is implemented 
using an auxilliary field Aj(r) which acts as a Lagrange multiplier. 



A key observation now is that the partition function in Eq. (6.10) is invariant under a 



site, i, and r-dependent U(l) gauge transformation Qir) where the fields transform as in 



Eq. (6.7), and A transforms as 



Xi — > Xi — 



dr 



(6.11) 



In other words, A transforms like the temporal component of a U(l) gauge field. 

How do we obtain the spatial components of the gauge field? For this, we apply the 



Hubbard-Stratonovich transformation of Eq. (2.23) to the tij hopping term in Eq. (6.9). 
For this, we introduce another auxiliary complex field Qijir) which lives on the links of the 
triangular lattice and replace the hopping term by 



E 



Q^j{r)fLfja-QUT)e 



(6.12) 



We now see from Eq. (6.7), that Qij transforms under the gauge transformation in Eq. (6.7) 



as 



Qij ^ Qij^^^ ^■'^ . 



(6.13) 



In other words, axg{Qij) is the needed spatial component of the compact U(l) gauge field. 
So far, we have apparently only succeeded in making our analysis of the Hubbard model in 



38 



Eq. (1.1) more complicated. Instead of the functional integral of the single complex fermion 



we now have a functional integral over the complex fermions fia, the rotor "i?,, and the 
auxilliary fields Aj and Qij. How can this be helpful? The point, of course, is that the new 
variables help us access new phases and critical points which were inaccessible using the 
electron operators, and these phases have strong correlations which are far removed from 
those of weakly interacting electrons. 

The utility of the new representation is predicated on the assumption that the fluctuations 
in the auxiliary fields Qij and Aj are small along certain directions in parameter space. So 
let us proceed with this assumption, and describe the structure of the phases so obtained. 
We parameterize 



Q^ 



A,; 



-iX — A- 



(6.14) 



and ignore fluctuations in the complex numbers Q^j, and the real number A. With these 



definitions, it is clear from Eqs. (6.11) and (6.13) that Aij and A^- form the spatial and 



temporal components of a U(l) gauge field, and so must enter into all physical quantities 
in a gauge invariant manner. The values of and A are determined by a suitable saddle- 



point analysis of the partition function, and ensure that the constraint (6.6) is obeyed. With 



these assumptions, the partition function separates into separate fermionic and rotor degrees 
of freedom interacting via their coupling to a common U(l) gauge field {Air, Aij). In the 
continuum limit, the gauge fields become a conventional U(l) gauge field A^ = {Ar,A). 
The partition function of the gauge theory is 

2 = 1 Vf,^{T)Vfl{T)VMT)Vn„{T)VAUr)VA,,{T) 
exp ^— J dr Cf + Cr + i A. 

^/ = E + ^ - ^^-) - E Q^^fl^^''f: 

i ^ ij 

£^ = -z ^ n^i (^-^ - Air^ - ^ QijC 



+ ^ l-Xriri - n{nri + 1) + 



u 



ririiUri + 1) I • 



(6.15) 



Thus we have fermions fia moving in a band structure which is roughly the same as that of 



the electrons in Eq. (6.1), rotors which obey a boson Hubbard-like Hamiltonian, and both 



are minimally coupled to a compact U(l) gauge field. 

We begin by neglecting the gauge fields, and computing the separate phase diagrams of 
Cf and Cr- 

The fermions are free, and so occupy the negative energy states determined by the chem- 
ical potential A. 
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FIG. 15: Possible phase diagram of the electron Hubbard model in Eq. (1.1) on the triangular 
lattice. This phase diagram is obtained by a mean-field analysis of the theory Cr in Eq. (6.15), 
similar to that for the boson Hubbard model in Chapter 9 of Ref. [19l Only the Mott insulating 
lobes with = —1,0,1 are compatible with the constraint in Eq. (6.6); these Mott insulating 

= + 1. 



lobes have fermion density ( 



The phase diagram of is more interesting: it involves strong interactions between the 
rotors. It can be a analyzed in a manner similar to that of the boson Hubbard model (see 



Chapter 9 of Ref. [T9l) . leading to the familiar "Mott lobe" structure shown in Fig. 15 



At large values of Q/f/ we have the analog of the superfluid states of the boson Hubbard 
model, in which there is a condensate of the rotor ladder operator e*'^. However, this op- 
erator is charged under the U(l) gauge field, and so this phase does not break any global 
symmetries. Instead it is a Higgs phase, like the Neel phase in the model of Section IV In 
the presence of the Higgs condensate, the operator relation in Eq. (6.4) implies that ~ /«, 
and so the fa fermions carry the same quantum numbers as the physical electron. Conse- 
quently, the fa Fermi surface is simply an electron Fermi surface. Furthermore, the Higgs 
condensate quenches the fluctuations, and so there are no singular interactions between 
the Fermi surface excitations. This identifies the present phase as the familiar Fermi liquid, 



as identified in Fig. 15 



Having reproduced a previously known phase of the Hubbard model in the U(l) gauge 
theory, let us now examine the new phases within the 'Mott lobes' of Fig. 15 In these states, 
the rotor excitations are gapped, and the rotor angular momentum has integer expectation 
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values. The constraint in Eq. (6.6) implies that only n,. = —1,0,1 are acceptable values. 



and so only these values are shown. It is clear from the representation in Eq. (6.5) that any 



excitation involving change in electron number must involve a rotor excitation, and so the 
rotor gap implies a gap in excitations carrying non-zero electron number. This identifies 
the present phases as insulators. Thus the phase boundary out of the lobes in Fig. 15 is a 
metal-insulator transition. 

Section |VIIB| will present an explicit demonstration of the insulating and metallic prop- 



erties of the phases in Fig. [15] by a computation of the transport properties of a broader 
class of models. 



The three insulators in Fig. 15 have very different physical characteristics. 



Using the constraint in Eq. (6.6) we see that the rir = —1 insulator has no fa fermions. 



Consequently this is just the trivial empty state of the Hubbard model, with no electrons. 

Similarly, we see that the rir = 1 insulator has 2 fa fermions on each site. This is the 
just the fully-filled state of the Hubbard model, with all electronic states occupied. It is a 
band insulator. 

Finally, we turn to the most interesting insulator with = 0. Now the electronic states 
are half-filled, with (flfa) = 1- Thus there is an unpaired fermion on each site, and its 
spin is free to fiuctuate. There is a non-trivial wavefunction in the spin sector, realizing an 
insulator which is a 'spin liquid'. In our present mean field theory, the spin wavefunction is 
specified by Fermi surface state of the /„ fermions. Going beyond mean-field theory, we have 
to consider the fiuctuations of the gauge field, and determine if they destabilize the spin 



liquid, as we did in earlier Section IV Here the fa fermions carry the A^, gauge charge, and 



these fermions form a Fermi surface. This is a crucial difference from Section [rV] where the 
iIj± fermions were gapped. In Section IV we found that the monopoles proliferated, leading 



to confinement and VBS order. Here, the gapless fermionic excitations at the Fermi surface 
prevent the proliferation of monopoles: the low energy fermions suppress the tunneling event 
associated with global change in fiuxP'^. Thus the emergent U(l) gauge field remains 
in a deconfined phase, and this spin liquid state is stable. These gapless gauge excitations 
have strong interactions with the fa fermions, and this leads to strong critical damping of 
the fermions at the Fermi surface which is described by a strongly-coupled field theorj^^^^. 
The effect of the gauge fluctuations is also often expressed in terms of an improved trial 
wavefunction for the spin liquid^-: we take the free fermion state of the fa fermions, and 
apply a projection operator which removes all components which violate the constraint in 



Eq. (6.6). This yields the 'Gutzwiller projected' state 



I spin liquid) = j Y\ 



1 — ( — J^^Ea/ia/" 



n ) io), 

\k<kp 



(6.16) 



where the product over k is over all points inside the Fermi surface. 



Finally, we turn to an interesting quantum phase transition in Fig. 15 This is the 
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transition between the spin liquid and the Fermi hquid at total electron density N" = 1, 
which occurs at the tip of the = Mott lobe. From the rotor sector, this looks like a 
Higgs transition, of the condensation of a complex scalar in the presence of a fluctuating 
U(l) gauge field. However, the fermionic sector is crucial in determining the nature of this 
transition. Indeed, in the absence of the Fermi surface, this transition would not even exist 
beyond mean field theory: this is because the U(l) gauge field is compact, and the scalar 
carries unit charge, and so the confining and Higgs phases of this gauge theory are smoothly 
connected. So we have to combine the Higgs theory of a complex scalar with the gapless 
Fermi surface excitations. Introducing a Bose field 

b = e-'^, (6.17) 

and coarse-graining to a continuum limit for the bosons, we find the field theorj^ 



C = \id^ + iA^)b\^ + s\b\^ + u\b\^ + tAM 



+ fi 



fa 



(6.18) 



where the energy e/ is to be adjusted to yield total fermion density M = 1. The transition 
is accessed by tuning s, and we move from a spin liquid for s < Sc, to a Fermi liquid for 
s > Sc- The critical properties of the theory in Eq. (6.18) have been studiecP^'^, and an 



interesting result is obtained: the Fermi surface excitations damp the gauge bosons so that 
they become ineffective in coupling to the critical b fluctuations. Consequently, the gauge 
bosons can be ignored in the b fluctuations, and the transition is in the universality class of 
the 2+1 dimensional XY model. 



VII. FRACTIONALIZED FERMI LIQUIDS 



In Section VI we met the canonical description of a compressible metallic state, the Fermi 
liquid. The is the state adiabatically connected to the metallic state of non-interacting 
electrons. It has long-lived fermionic quasiparticle excitations along the Fermi surface, and 
the area enclosed by this Fermi surfaces obeys the Luttinger theorem. 

Here we explore an extended model which allows for other compressible phases of strongly 
interacting electrons at generic densities which do not break any global symmetries, and 
which are not adiabatically connected to the limit of non-interacting electrons. We shall 
focus here on the compressible stat^^El i^^own as the fractionalized Fermi liquid (FL*). In 
principle, the FL* state can appear in a variety of models of correlated electrons, including 
ones with a single band, and all the sites equivalent with Ui = U. Such single-band FL* states 
have been described in recent workP^. However, these single-band analyses are involved, 
and require intermediate steps which make them sub-optimal for a first description of the 
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FIG. 16: The bilayer triangular lattice. The top layer (A) has nearest neighbor hopping t^, the 
bottom layer (B) and nearest neighbor hopping ts, and the inter-layer hopping is w. A closely 
related model is realized in the experiments of Ref. [6l 

FL* state. 

Instead, we will introduce the FL* state in a model with 2 types of inequivalent sites. 
As a simple example, consider the Hubbard model on a bilayer triangular lattice shown in 
Fig. [16] We label the two layers as A and B, and so there are 2 electron operators, CAia and 
CBia- We write the Hamiltonian as 

H = Ha + Hb + Hab 

Ha = -tA ^ C^iaCAja + H.C. + (c^ - 1^)"^ ('^Ait + ^AiJ,) 



(ij) 



Hp 



Hab 



+ UaY^ (^nAi^ - - J yiAti 

-tB ^ C^Bia'^Bja + H.C. + (e^ - /i) ^ (^Bit + ^^^4) 

+ f^B 5Z ( ^^^"^ ~ 
i ^ 

^AiaCBia + H.C. 



nBii 



(7.1) 



Here the sites i,j lie on a triangular lattice, and (ij) represents the sum over nearest- 
neighbor pairs. The Hubbard models on the two layers have distinct values of the hopping 
parameters, on-site repulsion, and on-site energies eA,B- Finally, there is an on-site inter- 
layer tunneling, w. Experiment^ on bilayer films of ^He adsorbed on graphite provide a 
remarkable realization of a closely related model. 

First, let us discuss the FL state, where Ua,b can be treated perturbatively. Diagonalizing 
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the one-electron Hamiltonian, we find two bands corresponding to the bonding and anti- 
bonding states between the two layers. Let be the total density of electrons for each 
bilayer site of the triangular lattice. So 

E {{^AaCAa) + (cLcBa)) = M. (7.2) 
a 

This relation holds for every site i, and the site-index has been left implicit. Depending 
upon the value of TV and interlayer tunneling one or both of the bands will be occupied, 
leading to one or two Fermi surfaces. Let the areas enclosed by the Fermi surfaces be Ai and 
A2'i if there is only one Fermi surface, A2 = 0. Luttinger's theorem fixes the areas of Fermi 
surfaces to a value which is independent of the nature of the electron-electron interactions. 
There is one Luttinger theorem for each global U(l) symmetry of the Hamiltonian which is 
not spontaneously broken in the ground stat d^ * ^^ *^. Here, the total numbers of both up-spin 
and down-spin electrons are separately conserved, and so there are 2 Luttinger constraints. 
However, we will implicitly only consider states in which spin rotation invariance is preserved, 
and so there is only a single constraint. The constraint is the same as that for non-interacting 



electrons, which, as in Eq. (6.3), is 

A1 + A2 



27r2 



M. (7.3) 



We will implicitly assume M > 1 below. 

We now wish to induce a quantum phase transition to a FL* state. This is most easily 
done in a model in which the bottom layer B has a density of one electron per site, while 
the top layer A remains dilute (as in the experiment in Ref. |6]). For small hopping this is 
achieved ioi Ua = Ub = U and 

eb <eA<eB + U . (7.4) 



Then the bottom layer B will acquire strong electronic correlations like those in Section VI 



while the dilute gas on layer A can be treated perturbatively in the two-particle scattering 



amplitude. It is customary at this point to follow the analysis of Section |II C[ and project 
onto this restricted Hilbert space, while using a canonical transformation to derive an ef- 
fective Hamiltonian. The restricted space has only spin degrees of freedom on B lattice 



sites, and as in Section II C, these spins have exchange interactions with each other. The 
canonical transformation also generates exchange interactions between electrons separate 
layers, and this is known as the Kondo exchange interaction. The resulting Hamiltonian 
is the Kondo- Heisenberg model. However, we will not take this step here, and continue to 



work with the Hubbard model in Eq. (7.1). 



We will assume that layer 5, with a density of one electron per site, realizes the spin 



liquid state discussed in Section VI ; i. e. it is the = spin liquid in Fig. 15 , with a spinon 



Fermi surface. We can obtain a description of this spin liquid by applying the analysis of 
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Section VI to layer. So we replace Eq. (6.4) by 



csa = e 



(7.5) 



and perform the same transformations which led to Eq. (6.15). Then we take the same 



continuum limit as that used for Eq. (6.18), and obtain the following continuum Lagrangian 



L which captures the low energy physics of the Hubbard model in Eq. (7.1 ). The degrees of 



freedom are the A layer electrons c^a, the B layer spinous /„, and the bosonic rotors b = e 



as in Eq. (6.17). The structure of the terms also follows from general considerations of gauge 



invariance and the preservation of global symmetries. 



d . . 1 



. tf -iAr- 
OT 2m f 



fa 



{df, - {tr - fi)S^r - iA^ + iAcxt,fi) 

{d^ + {tr - fi)S^r + iA^ - iA,^t,^^) b] + s\b\^ + u\b\^ 



^Aa 



d 1 



(7.6) 



Here A^ = {A^-, A) is an emergent U(l) gauge field; we have also introduced a non-fluctuating 
electromagnetic gauge fleld A(,xt,^l as a source term which couples to the current of the 
globally conserved electromagnetic charge; we have coarse-grained 6 to a complex scalar 
field with both amplitude and phase fluctuations; the symbol /i refers separately to the 
chemical potential and spacetime component, and the interpretation should be clear from 
the context; the flnal Yukawa term is allowed by the symmetries, and represents the inter- 
layer tunneling w; the on-site energies e/ and er are related to and and have to be 
tuned so that the system obeys the density constraints to be discussed below. 



To review, the continuum theory in Eq. (7.6) has a U(l) xU(l)cxt symmetry associated 
with the transformations 



fa ^ /aC 
fa fa 



b ^ b , CAa^ CAa 

b b'^ , CAa CAaC'^ 



(7.7) 



The flrst U(l) symmetry is gauged by the dynamical emergent U(l) gauge fleld A^, and is 



the same as that in Eq. (6.7). The second U(l) symmetry remains global; the flxed external 
electromagnetic fleld Acxt,ii couples a source term which gauges this global symmetry. 

In general, there will be 2 Luttinger constraints associated with these two U(l) 
symmetrie d^ * ^^ * ^^ ! (as before, we are ignoring spin rotation symmetries here, which is as- 
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sumed to be always fully preserved). The first transformation in Eq. (7.7) leads to a Lut- 
tinger constraint on the associated conserved charge density (which is the continuum analog 



of Eq. (6.6)) 



27r2 



^fB. 



{7.1 



Here N'b is the density of electrons on layer B in the projected Hilbert space: our present 



lattice derivation was for A/b = 1, but the continuum theory in Eq. (7.6) is sensible for any 



value of N'b- The operator Qi, is the rotor angular momentum, given by 



Qb 



dfj. 



(7.9) 



Thus there must be a Fermi surface enclosing area Ai , which counts the density of / fermions 
minus the bosonic rotor density. 



Similarly, the second transformation of Eq. (6.6) leads to the constraint 



U-Mb. 



(7.10) 



Again, there is a Fermi surface enclosing area A2 which counts the density of ca fermions 
minus the bosonic rotor density. Thus our analysis so far appears to imply that there must 
be at least 2 Fermi surfaces, and their areas are constrained by the two independent relations 



in Eqs. (7.8) and (7.10) 



This last conclusion seems rather surprising from our discussion above of the FL phase. 



There we found only a single constraint in Eq. (7.3) for the total areas of one or more 



Fermi surfaces. The only possible conclusion is that the FL phase is not one in which the 



U(l) xU(l)ext symmetry of the Lagrangian C in Eq. (7.6) remains unbroken. Rather the 
FL phase is realized as a Higgs phase in which the U(l)xU(l)cxt symmetry in Eq. (7.7) is 



broken down to a diagonal U(l). Just as in Section |VI| this is the Higgs phase in which the 
boson h condenses 

(6) 7^ in the FL phase. (7.11) 

Once the symmetry is broken in this manner, the corresponding Luttinger constraint no 
longer applied^'^^'^. Only the sum of the constraints in Eqs. (7.8) and (7.10) applies, and 



this leads immediately to the defining relation in Eq. (7.3 ) of the FL phase. The condensation 



of h also quenches the emergent U(l) gauge field, so there are no gapless gauge excitations 



in the FL state, again as in Section VI 



We now see that state of the theory C in which the Luttinger constraints in Eqs. (7.8) 



and (7.10) apply separately is a new phase: this is the advertised FL* phase, in which the 

(6) = in the FL* phase. (7.12) 



boson h is uncondensecP^ 
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The full U(l)xU(l)ext symmetry is preserved, and the gauge boson becomes an emergent 
gapless photon. The arguments for the stability of the FL* phase towards gauge fluctuations 
mirror those of Section |VI| for the stability of the spinon Fermi surface in the spin liquid. 



The criteria in Eqs. (7.11) and (7.12) show that the transition between the FL and 
FL* states is tuned by varying the coupling s in Cb from negative to positive values. The 
transition between these phases occurs at a quantum critical point where the scalar b is also 
critical. 



A. Connections to holographic metals 

We now connect the above generic theory of the compressible FL and FL* phases of 
the Hubbard model to recent studies of compressible metallic phases via the AdS/CFT 
correspondence. The discussion below refers to recent work from the gravity perspective; an 
analysis starting the canonical supersymmetric gauge theories of gauge-gravity duahty may 
be found in Ref. [51 

A connection was made in Ref. H] bet ween a mean-field solution of models like the Hubbard 



model in Eq. ( 7. 1 ) and a particular AdS realization of a holographic metal. Specifically, 



the bilayer Hubbard model has been solved in a limit with infinite-range hopping matrix 



elements between the sites (in contrast to the nearest-neighbor hopping shown in Eq. (7.1 )). 
A detailed correspondence was found between the low energy properties of the FL* phase 
of such a model and the holographic theor}J21IZlE5l which the low energy limit factorized 
to a AdS2xR'^ geometry [d is the dimensionality of space). This work has been recently 
reviewed in the companion articl^. 



However, the mean- field solution of Eq. (7.1) and the AdS2xR geometry share a number 



of artifacts: they have a non-zero ground state entropy, and the spin correlations of layer B 
scale with dynamic exponent z = oo. These artifacts are not expected to be properties of 



the field theory C in Eq. (7.6), applicable for models with short-range interactions. 

It is clearly of interest to move beyond the AdS2xR'^ factorization in the holographic 
theory, and derive a holographic model which has a closer correspondence with the phases 
of the field theory in Eq. (7.6). A number of recent theorie^^^'^ have examined the feedback 



of the finite density matter on the metric of the AdS space, and found that the AdS2 horizon 
disappears at T = 0, and is replaced by a metric with a finite value of z. Many physical 



properties of such holographic metals are similar to those of the field theory in Eq. (7.6), 
but a detailed correspondence awaits future work. 

It is useful to consider these recent works in the context of a holographic R( j74 | 75 | 8 5rilE] 
these works, the UV degrees of freedom are coupled to external sources, which in our case 
is v4cxt,^- Then an effective action is derived which couples the external sources to the IR 
degrees of freedom. Two distinct fixed-point theories have been considered in the literature: 
those of Nickel and Son!^, and of Faulkner et al^ and Faulkner and PolchinskP^. We argue 



47 



here that these fixed points capture the physics of the FL and FL* phases respectively. 

Let us consider, first, the theory of Nickel and Sori^. They argued that the low energy 
theory had an emergent U(l) gauge field, so that the theory had a U(l)giobaixU(l)gauge 
symmetry. This is strikingly similar to the U(l) xU(l)cxt symmetry of Eq. ( |7.6 ). Indeed, 



we can more closely map the low energy theory of the FL phase of Eq. (|7.6|) to the model 
proposed by Nickel and S( 
fiuctuations of its phase b = e~ 



proposed by Nickel and Son. In the FL phase, we condense the b boson, and focus on the 

xhen the effective theory of the FL phase of Eq. (7.6) is 



(7.13) 



where IIj is the effective action obtain after integrating out the / spinous. The structure of 
Eq. (7.13) is essentially identical to Eqs. (6) and (52) of Nickel and Sonl^Sl. 

Consider, next, the corresponding low-energy theory of FL* phase. Now the b field is not 



condensed, and has an energy gap, A. So we can safely integrate it out from Eq. (7.6), and 
obtain an effective theory for the electrons, the spinous, and the gauge fields: 

+ [V{Ar - Aext,r) " 9,(A - A^xt)]' + ^^4 [V X (A - Aext)]' , (7.14) 

where Ji ~ J2 ~ w^/A. The coupling Ji is the Kondo exchange between the electrons in 
layer A and the spins on layer 5, while J2 couples density fiuctuations of the two layers. A 
key property of the FL* phase is that the couplings Ji,2 can be treated perturbatively: there 
is no fiow to strong coupling in the Kondo exchange, and the layer B spins are not screened 



by the conduction electrons. Let us now rewrite the matter component of Eq. (7.14) as 



FL* 



Fl_CAa + (^Aa-Fa 



(7.15) 



where Fa is a IR fermion defined bj3 



Fa ^ -Ji (Kpf^a'^^sfs) CAP - J2 {f%) CAc. 



(7.16) 



Notice that both fermions in the displayed term in Eq. (7.15) are invariant under the emer- 



gent U(l); this term is a coupling between the microscopic fermion CAa and a composite 
gauge-invariant fermion operator F^ representing the IR degrees of freedom. We can view 



Fa in Eq. (7.16) as the most general fermion operator which involves the IR fermions /„, 



which is invariant under the gauge transformation associated with A^, and which also carries 



the global electron number charge associated with Aext,^- Then structure of Eq. (7.15) is pre- 
cisely that of the semi-holographic theory of Faulkner et al^ and Faulkner and PolchinskP^l^ 
and their IR fermion Fa is chosen by essentially identical criteria. 

It would clearly be of interest to also find another fixed point of the holographic theory 
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corresponding to the quantum-critical points between the FL and FL* phases. 



B. Transport theory 



We conclude our discussion of FL and FL* phases by presenting a general formulation 
of their transport properties. The arguments below are in the spirit of those of loffe and 
LarkiiP. 



We begin with a theory like £ in Eq. (7.6), and integrate out the matter fields to obtain 
a Coleman- Weinberg effective action for the U(l)xU(l) gauge fields and Aext,/^- In 
general, the form of this effective action is constrained only spatial isotropy and gauge 



invariance. Using the projectors defined in Eq. (3.13), we can write the quadratic portion 



of the effective action in the following form (we work here in Euclidean time, and a;„ is a 
Matsubara frequency) 



S 



+ {A^ - Axt,^) {P;:,Kj;{uJn, k) + Pl,K^{uJn, fc)) (^M ' ^e.t,^) 

+ A.t,;.(Pii^e^(c^n, k) + PlKl{uj^, A:)) A.t,, 



(7.17) 



Here i^f are given by the correlator of the current of the / fermions, K^'"^ by the correlator 



L T 

of the current of the ca fermions, and K^^ ' by the current of the bosonic b rotors. Note 



in the definition of the 



that, unlike Eq. (3.12), we have not pulled out a factor of \Jk'^ + 
K^'^ . In general, determining these functions requires a complex transport analysis of the 



theory in Eq. (7.6). However, in the FL and FL* phases, the simpler low energy effective 
theories in Eqs. (7.13) and (7.14) lead to simple forms for the bosonic correlators K, ' . 



In the FL phase, integrating out the phase in Eq. (7.13) we obtain 



KU^n,k) 



Kl{L0n,k)=K2. 



(7.18) 



Thus K, ' are constants in the limits of small momenta or frequency. Indeed, had we chosen 



the velocity of 'light' judiciously in the definition of P^^ in Eq. (3.13), we would have had 



In contrast, in the FL* phase, we can directly match the low energy theory in Eq. (7.14) 



to Eq. (7.17) and obtain 



Kt{un,k) = Ks{e + ul 



K^iun,k)=Ksul + K,k' 



(7.19) 
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Now the K^^ ' vanish in the hmit of small momentum and frequency. 

We need to use the respective low energy theories of the FL and FL* phases in Eqs. (7.13) 



and ( 7.14[ ) to determine Kj^'^ and K^''^, and then combine them with the above results for 



-K"^'"^ to obtain the physical conductivity. As in Nickel and SoeP^, and in loffe and Larkinpl, 
this is obtained by implementing the equation of motion of the emergent gauge field Afj_ in 



Eq. (7.17). This equation of motion is equivalent to the constraint that the current of the 



b bosons must equal the current of the / fermions, which is a consequence of the lattice 



constraint in Eq. (6.6). Evaluating the equation of motion 5S/SA^ = from Eq. (7.17) 



and substituting the resulting value of back (after suitable gauge fixing), we obtain an 
effective action for the probe field Aext,/^ alone 



ext 



k) + PlKl,{uj^,k)]A^ 



(7.20) 



with 



Klti^n,k) = Kj{un,k) + 



K}{un,k)Kl;{un,k) 
Kf{un,k) + Kl;:{un,k) 

Kj{ujn,k)K^{Un,k) 
Kj{un,k) + Kf{ujn,k) 



(7.21) 



After analytic continuation to Minkowski space, these results lead to the physical conduc- 



tivity via the Kubo formula in Eq. (2.33) 



00 



(7.22) 



The distinction from Eq. (3.15) is a consequence of omitting in Eq. (7.17) the prefactor 



iyaJ|~+~P present in Eq. (3.12). 

These expressions can be used along with specific computations of the dynamics of the 
/ and c fermions: the latter can be carried out either using a Boltzmann theory of the 
continuum model in Eq. (7.6), or via a theory on AdS. The analysis by Nickel and SorP^for 



their holographic theory is equivalent to the application of Eq. (7.21). Let us verify that the 



present method yields the expected FL behavior in the Higgs phase where Eq. (7.18) implies 
that K^{u,0) = K2- We assume that c and / Fermi surfaces have metallic conduction with 
K^{uj,Q) = —iucTc and Kj{u,0) = —iuaf, with acj the respective conductivities. Inserting 



these expressions in Eq. (7.22), we obtain the expected FL behavior with a = ac + Cf in the 



limit a; — )■ 0. Thus there is no superfluidity associated with the condensation of b, and the 
gauge fluctuations lead eventually to metallic behavior. Similarly, it is easy to show that 
if b and c excitations are gapped, we have insulating transport, even though the / spinous 
have a gapless Fermi surface. 

Further theoretical work exploring the connection between the AdS and Boltzmann ap- 
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proaches to transport is clearly of interest. 
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